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^ . Abstract 

We develop further the theory of Rational Conformal Field Theories (RCFTs) on a 
cylinder with specified boundary conditions emphasizing the role of a triplet of alge- 
bras: the Verlinde, graph fusion and Pasquier algebras. We show that solving Cardy's 
equation, expressing consistency of a RCFT on a cylinder, is equivalent to finding inte- 
ger valued matrix representations of the Verlinde algebra. These matrices allow us to 
naturally associate a graph G to each RCFT such that the conformal boundary condi- 
tions are labelled by the nodes of G. This approach is carried to completion for sl{2) 
theories leading to complete sets of conformal boundary conditions, their associated 
cylinder partition functions and the A-D-E classification. We also review the current 
status for WZW sl(3) theories. Finally a systematic generalization of the formalism 
of Cardy-Lewellen is developed to allow for multiplicities arising from more general 
representations of the Verlinde algebra. We obtain information on the bulk-boundary 
coefficients and reproduce the relevant algebraic structures from the sewing constraints. 
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1 . Introduction 



1.1 History and motivation 



The subject of boundary conformal field theory has a fairly long history. It was born more 
than ten years ago, in parallel work on open string theory |TJ, ^ [| f|, |J| [], [7], |S| and on 
conformal field theories (CFTs) describing critical systems with boundaries ||. The work 
of Cardy |10| was a landmark, leading to the unification of methods, to the introduction of 
important concepts such as boundary conformal fields and to the systematic investigation of 
their properties and couplings [p| , [12]] . The subject remained dormant for some time, in spite 



of some activity motivated again by string theory |13| [14], [L5| and of beautiful applications 
to the Kondo problem |16|, [T7|. Lately, the subject has undergone a revival of interest in con- 
nection with various problems. On the one hand, work on boundary conditions in integrable 
field theories and boundary flows |Tj| [T9|, ^] and on quantum impurities [|2T], E2| motivated 
a closer look at boundary CFT. On the other hand, within statistical mechanics, integrable 
Boltzmann weights satisfying the so-called Boundary Yang-Baxter Equation (BYBE) were 
constructed in lattice models |23|, |24|. Finally new progress in string theory was another 
reason to reconsider the problem. Generalizations of D-branes as boundary conditions in 
CFT have been studied by several groups |^5j | 



27], H, ||, |0|. 



In the present work, we want to reconsider several issues in the discussion of boundary 
conditions in (rational) conformal field theories: what are the general boundary conditions 
that may be imposed, what are the structure constants of the bulk and boundary fields in 
the presence of these boundary conditions. These are the basic questions that we want to 
address. The methods that we use are not essentially new, but are based on the systematic 
exploitation of the work of Cardy and Lewellen [pi], [TT|, [12 . 

Among the main results of this paper: 

• We establish a connection between the classification of boundary conditions and the 
classification of integer valued representations of the fusion algebra. A preliminary 



account of this result was given in |31], |32[] . In the same vein, we show that it is natural 
to associate graphs to these problems. In particular, an ADE classification of boundary 
conditions for Wess-Zumino-Witten (WZW) and minimal s/(2) theories emerges in a 
natural and simple way. A discussion of the state of the art for sl(3) models is also 
included. 

We point out the deep connections between the features of conformal field theory in 
the bulk and in the presence of boundaries. The classification of the latter has some 
bearing on the classification of the bulk properties (modular invariants etc). This is 
not a new observation. In particular, in string theory many connections are known 
to exist between open and closed string sectors, but it seems that the point had not 
been stressed enough. A triplet of algebras, specifically the graph fusion algebra and its 
dual, the Pasquier algebra, appears naturally in our discussion, along with the Verlinde 
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algebra. 



We reanalyse in a systematic way the couplings (structure constants) of fields in the 
presence of boundaries and the equations they satisfy, generalising the formalism of 
Cardy-Lewellen to accomodate the appearance of nontrivial multiplicities. In the di- 
agonal cases we find a direct relation between the chiral duality identities of Moore 
and Seiberg and the basic sewing relations of the boundary CFT. The main point is 
the observation that the bulk boundary coefficients in the diagonal case essentially 
coincide with the matrices S(j) of modular transformations of torus 1-point functions. 



In this way the two basic bulk-boundary equations [11], |T2| are shown to be equivalent 
to the torus duality identity of |33| . 



Some more particular results include the extension of Cardy's equation to non-specialized 
characters, thus lifting an ambiguity in the original derivation, the proof of uniqueness of 
boundary conditions for sl(2) WZW and minimal models and sl(N)i models, the clarifica- 
tion of the role of the graph algebra and the recovery of this algebra along with its dual, the 
Pasquier algebra, from the boundary sewing constraints. 

1.2 Background on bulk CFT 

In this paper, we are only concerned with Rational Conformal Field Theories (RCFTs) . We 
first establish notations etc. In the study of a RCFT, one first specifies a chiral algebra 
A. It is the Virasoro algebra or one of its extensions: current algebra, W algebra etc. The 
generators of this algebra will be denoted generically W n and include the Virasoro generators 
L n . At a given level, the theory is rational, i.e. A has only a finite set X of irreducible 
representations Vj, i EX. The label i* indexes the representation conjugate to i, and i — 1 
refers to the identity (or vacuum) representation. We also suppose that the characters 
Xi(l) = tr Vj<? L °~^ of these representations, the matrix S of modular transformations of the 
x's and the fusion coefficients iVy of the V's are all known. The matrix 5y is symmetric 
and unitary and satisfies S 2 = C, where C is the conjugation matrix Cy = 5ji*. The fusion 



coefficients are assumed to be given in terms of S by the Verlinde formula [34 



an assumption that rules out some cases of RCFTs. 

A physical conformal theory is then defined by a collection of bulk and boundary fields 
and their 3-point couplings (OPE Coefficients). In particular, the spectrum of bulk fields 
is described by the finite set Spec of pairs (j, j) of representations, possibly appearing with 
some multiplicities Njj, of the left and right copies of A, such that the Hilbert space of the 
theory on an infinitely long cylinder reads 

W = eyj )6 s P ecV i ®Vj , (1.2) 
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with the same multiplicities Njj. The modular invariant torus partition function 

Ztorus = E N fj Xj(q) (xj(qj) (1-3) 



is a convenient way to encode this information. The finite subset £ of labels of elements of 
the spectrum that are left-right symmetric will play a central role in the following 

£ = {j\(j,3 = j) eSpec) , (1.4) 

and will be called the set of exponents of the theory. Recall that these exponents may come 
with some multiplicities. To distinguish them as different elements of the set £ a second 
index will be often added, i.e., (j, a) G £ , for j G X. 

In terms of all these data, one is in principle able to compute exactly all correlation 



functions of the CFT on an arbitrary 2D surface, with or without boundaries [35, [33], 12]. 
These data, however, are subject to consistency constraints: single-valuedness of n-point 
functions on the plane, modular invariance of the torus or annulus partition function, etc, 
all rooted in the locality properties of the theory. In this paper, we shall reexamine the 
conditions that stem from surfaces with boundaries (half-plane or disk, cylinder or annulus) 
and explore their consequences. 

For later reference, let us also recall that RCFTs fall in two classes. In the first class 
("type I"), the Hilbert space ( |1.2|) is a diagonal sum of representations of a larger, "extended" , 
algebra A' D A. Accordingly, the partition function ( |1.3| ) is a sum of squares of sums of 
characters 

z= E IEx.1 2 - (i-5) 

blocks B ieB 

The second class ( "type II" ) is obtained from the first by letting an automorphism ( of the 
fusion rules of the extended algebra A' act on the right components, thus resulting in a 
non-block-diagonal partition function 

^ = E(Ex,)( E xj). (i.6) 

B ieB je((B) 

For example, in the classical case of sl(2) theories, classified by A-D-E Dynkin diagrams, 
the A, D2 P , Eq and Eg cases are of the first type, whereas the D2 P +\ and E7 are obtained 
respectively from the A 4p _i and D w cases by a Z 2 automorphism of their fusion rules. We 
shall see below that the study of boundary conditions on a cylinder has some bearing on 
these expressions of torus partition functions. 



2 . Cardy Equation and Verlinde Algebra 
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2.1 Boundary states 



As discussed in fl36| , on the boundary of a domain such as the upper half plane or a semi- 
infinite cylinder, one must impose a continuity condition of the form 

T{z) = T{z)\ z=z W(z) = W(z)\ z=2 . (2.1) 

While the first of these conditions has the direct physical meaning of the absence of energy- 
momentum flow across the boundary, or the preservation of the real boundary by diffeo- 
morphisms, the condition(s) on the other W may be generalized to incorporate a possible 
"gluing automorphism" BT], |25|, |27| 



W{z)=nW{z)\ z=z . (2.2) 

A semi-annular domain in the upper half-plane may be conformally mapped into an annulus 
in the complex plane by C, = exp(— 2iirw/T), w = ^logz. Then as shown by Ishibashi [Q] 
and Cardy flOfl , the boundary condition becomes 

( 2 T{() = fT(C) ( SW W{() = (-() SW W(0 for |C| = 1 and |C| = e 2wL/T (2.3) 

where sw denotes the spin of W, or more generally 

C w w(() = (-{) sw nW{() . 

Through radial quantization, this translates into a condition on boundary states \a)^ 

(W n - (-l) Sw n(W_ n )) \a) n = . (2.4) 

This includes in particular the condition that 

(L n -Z_ n )|a) n = 0, (2.5) 

assuming that the automorphism Q keeps invariant the Virasoro generators. For the central 
charge operator we have (k — k)\a)n = 0. 

Solutions to this linear system are spanned by special states called Ishibashi states 0, 
labelled by the finite set £q = j = ^(j)) £ Spec}, where ui depends in particular on 
Q. To see this, let us consider first the simpler equation ( |2.5| ) in the case when A is the 
Virasoro algebra and Q is trivial. Then observe that one may solve ( |2.5| ) in each component 
of ( |1.2|) independently, as these spaces are invariant under the action of the two copies 
of A. Now we recall that any state A = J2 n ,n a n,n\j,n) ® \j, n) in Vj ® Vj is in one-to- 
one correspondence with a homomorphism Xa = J2n,n a n,n\j, n\ of Vj into Vj. This 
uses the scalar product in Vj. Since L_ n = L\ for that scalar product, ( |2.5|) implies that 
L u Xa = Xj{L n , i.e. that X^ intertwines the action of L n in the two representations Vj 
and Vj of the Virasoro algebra. As these two representations are irreducible, they must be 
equivalent, which by our convention on the labelling of representations, means that j = j. 
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Thus the only non- vanishing components of A in Ti are in diagonal products Vj (g) Vj and in 
each one, Xa is proportional to the projector Pj = J2 n \ji n )(ji n \- To fix the normalization 
we choose Xa = Pj and the corresponding Ishibashi state is denoted This completes 
the proof that there is an independent boundary state \j)) for each element of the set 
£ = {j\(j,J =j) e Spec}. 

The argument is a formal extension of the proof, based on the Schur lemma, of the 
existence and uniqueness of an invariant in the tensor product of finite dimensional rep- 
resentations. It extends to the odd spin sw case (fOj). We have to use the fact that 

= (— l) Sw " U^ 1 W- n U with respect to a bilinear (or hermitian) form where U is a uni- 
tary (or antiunitary) operator. One exploits the same definition of the homomorphism 
Xa '■ Vj' — > Vj, now Vj> ,Vj being highest weight representations of the chiral algebra A 
generated by W n . However Xa corresponds to states in Vj ®Uq_U Vj', where Un is a unitary 
(antiunitary) operator implementing the automorphism Q(W n ) = Uq W n U^ 1 . The equation 
Q2.4Q leads to W n X A = XAW n again with the result f = j and Xa = Pj while the Ishibashi 
states are given by \j)) a = J2n \j, n) ®U n U \j, n). 

The operator U is in general non-trivial, e.g., for the sl(N)k WZW theories = 
w(W- n ) where w is the horizontal projection of the Chevalley involution of the affine algebra 
P8|] , i.e., it is determined by w(e ai ) = —f ai , w(f ai ) = —e ai , for the simple roots aij of 
sl(N), where e a / f a are raising/lowering operators respectively, and for the Cartan generators 
w{h l ) = —h % ,i = 1, 2, . . . , iV — 1. The pair (j, uj{j)) characterising £q refers to the eigenvalues 
of h l on the first term in \j))n- If ^ is the identity, then u(j) = —j. U Q = wq, where 
wo represents the longest element of the Weyl group, then u>o(j) = — j* and hence u>(j) = 
Wow(j) = j*. On the other hand u)(j) = j for Q coinciding with the Chevalley automorphism 
Q = w. In the last two of these examples we can identify Vj = U^U Vj with a highest weight 
module Vu(j) , w(J) £ X. It should be stressed that all these automorphisms fl keep invariant 
the Sugawara Virasoro generators so the condition ( |2.5| ) is also satisfied on the corresponding 
Ishibashi states. 

We shall hereafter drop the explicit dependence on Q. 

In fact we still have to define a norm (or a scalar product) on boundary states, in 
particular on the Ishibashi states. We have to face two difficulties. First because of the 
infinite dimension of the representation Vj, the most naive norm, proportional to tr Pj, 
would be infinite. The second problem concerns non-unitary representations. In such cases, 
the hermitian form on Vj is not positive definite, and we may encounter signs in the norm 
of states. 

The first problem requires some regularization of the naive norm. Let (p = e~ be 
a real number, < q < 1. Then ((j\q^ Lo+ ^°~^\j)) = tr Pjq L °~^ = Xj{o)- We write in 
general, allowing some multiplicity a = 1, ■ ■ ■ , Njj for the representations: 

((j',a'\q-> (Lo +L °-^\j, «» = %'<W XM) ■ (2-6) 

2 Many thanks to G. Watts (private communication) to whom we owe this elegant derivation. Some 
elements had appeared already in M. Bauer's PhD thesis (1989). 
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The norm of should then be some renormalized version of the q — > 1 limit of 
|fj. i.e., of the limit in which q = e 2nlT , the modular transform of q, tends to 0. In unitary 
theories, a new scalar product on boundary states may be defined according to 

({M\j'a')) = hmq^((f,a'\q^ Lo+L °-^\j,a)) 

= <%' 5 aa i Yaaq^xM) = $jj> <W<Sij , ( 2 -7) 

where we have used the fact that in a unitary theory, the leading character in the q — > 
limit is that of the identity operator Xi(q) ~ q~^ . Note that Sji is, up to a factor l/Sii, the 
quantum dimension of the representation j, a positive number. Thus in unitary theories the 
normalization chosen for Xj is such that the states \j)) are orthogonal for the scalar product 
( |2.7| ), with a square norm equal to S\j. Although in non- unitary theories the limit q — > 
in ( p.7|) does not exist in general, due to the existence of representations of conformal weight 
hi < that will dominate that limit, we may still define the norm by the same formula as 



2~7j ). Alternatively, if jo denotes the unique representation of smallest conformal weight 



hj < belonging to £, and c e fr := c — 24/i JO , then we may define 

«ja||jV» = Sjj.S^^Kmq^xM. (2.8) 

^303 q 

In all cases we thus have 

(0'«lbV}> = SjfSaa'Sxj (2.9) 

which is now of indefinite sign. In the sequel, we use more compact notations and the 
multiplicity label a will be implicit when referring to j £ £. 

The most general boundary state \a) satisfying condition fl2.4Q must be a linear combi- 
nation of these Ishibashi states, which, for later convenience, we write as 

l Q > = E-7= b» • (2-10) 

jes y oij 

We denote by V = {a} the set labelling the boundary states. We assume that an involution 
a — > a* in the set V is defined and that i\P a * = ip J a * = where j — > j* is an involution in 

£ (in general (j, a) — > (j*,a*), see Appendix B for examples). We define conjugate states 
as 

<&h=£«il-7=- (2-n) 



As explained in 25 , this may be regarded as resulting from the action of an 



antilinear CPT operation. As a consequence 



<%> = E ^FMOTO = MY (2-i2) 



jes Sl 3 jes 
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so that the orthonormality of the boundary states is equivalent to that of the if} 1 8. 

In some cases, such as in the computation of partition functions involving the specialised 
characters in the next section, it is sufficient to impose only the Virasoro condition (|2.5|) 
on the boundary states. Then the sum in (|2.10|) , when interpreted in terms of Ishibashi 
states pertaining to some extended symmetry, may include states \j))a with different Q. 
For example, in the minimal sl(2) models when multiplicities occur in S , one can build 
the Ishibashi states using the Coulomb gas realisation with A = u(l). Then there are two 
choices of Q keeping L n invariant. This allows, in particular, the construction of two different 
Ishibashi states with the same value of the scaling dimension, i.e., the explicit resolution of 
the degeneracy of states denoted \j, a)) above. Such mixtures of Ishibashi states may be 
used in determining the boundary states of the non-diagonal (A, D even ) models. 

2.2 The Cardy equation 

We now consider a conformal field theory on a finite cylinder. Following Cardy, the partition 
function may be expressed in two alternative ways. Regarded as resulting from the evolution 
of the system between boundary states a and b under the action of the Hamiltonian on (i.e. 
the translation operator along) the cylinder, it is 

Z b]a =(b\q3 {Lo+L °-* i) \a) (2.13) 

where <p = e~ 2n ^ describes the aspect-ratio of the cylinder of period T and length L as 

b 









A 

if 

\; , 
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\) 




a j 
(a) 




(b) 





Figure 1: The two computations of the partition function Zu a : (a) on the cylinder, between the 
boundary states a and b, (b) as a periodic time evolution on the strip, with boundary conditions a 
and b. 

in Figure 0. Decomposing the boundary states on the Ishibashi basis and using ( |2.6| ), one 
obtains 

£6[. = Ev£(Wr^ (2.14) 

where the states \j)) are admissible Ishibashi states of the system, i.e., the label j runs over 
the set £. 
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On the other hand, Zb\ a may be regarded as resulting from the periodic "time" evolution 
under the action of the translation operator along the finite width strip in the presence 
of boundary conditions a and h. The latter manifest themselves only in the nature of the 
Hilbert space TCb a and its decomposition into representations of a single chiral algebra: 7ib a = 

T 

©i^i6 a Vj with non-negative integer multiplicities riib a . If q = e _7r i, Zu a is a linear form in 
the characters 

Zb\a = J2^(q)n lb a ■ (2.15) 

We choose to write the modular transformation of characters in the form Xiio) = 
J2j SijXjio), hence Xj(<L) = 12i Sji'Xiiv) ■ Provided that specialized characters Xi{o) are 
considered, this complex conjugation is immaterial, since Xi(o) = Xi*(o)- We shall, however, 
make later use of unspecialized characters (Appendix A), for which it does matter. With 
this convention, and assuming for the time being the independence of characters, the two 



expressions (|2.14j ) and (|2.15|) are consistent provided 

"fa 6 = £ IK t# (W)* ■ (2- 16 ) 
jes ^3 

In the sequel we will refer to this as the Cardy equation. In the left hand side, we have used 
the first of the following symmetries 

n ia b = = n ib * a * , (2.17) 

which follow from the properties of the modular matrix and of the coefficients t/>£. 

The boundary states \a), \b), are thus such that n ia b is a non- negative integer. Unique- 
ness of the vacuum implies n la b < 1. The Cardy equation (2.16|) is a non-linear constraint 



on the components of boundary states \a) and \b) on the basis of Ishibashi states. Note 
also that it implies that J2i nib a Sji vanishes if j ; ^ £ and, except in cases with multiplicities, 
must factorize into a product of contributions of the a and b boundary states, a non-trivial 
constraint. Still, these constraints seem difficult to solve in full generality. 

Before we proceed, we have to pause on the question of independence of characters. 
In general, it is not true that specialized characters such as those that we have been using 
so far, are linearly independent. For instance, complex conjugate representations % and i* 
give rise to the same character Xi(<l)- Unfortunately, little is known about unspecialized 
characters for general chiral algebras, beside the case of affine algebras. In Appendix A, we 
show that in the case of rational conformal field theories with a current algebra, the previous 
discussion may indeed be repeated if the energy momentum tensor of the theory has been 
modified in such a way that unspecialized characters appear. Then using the known modular 
transformations of the latter ||38|| , one derives ( 12.161) . We shall therefore assume that (|2.16|) 



holds true for general RCFT. 

We now return to the Cardy equation fl2.16j ), and assume that we have found an or- 
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thonormal set of boundary states, i.e. satisfying 



E^(W)* = $*- (2-18) 

Moreover we make the stronger assumption that we have found a complete set of such states, 
i.e. satisfying 

E^'r = <w- (2.i9) 

a 

(Note that this implies that the number of these boundary states must be equal to the 
cardinality of £). 

Finally we recall that the ratios Sij/Sij, for a fixed j G X, form a one-dimensional 
representation of the fusion algebra, as a consequence of the Verlinde formula ( |1 . 1|) : 

EW^- (2-20) 



It follows from (|2.19|) , (|2.20| ) that the matrices Hi, defined by 

(n t ) a b = nj % E 1 (2.21) 



also satisfy the (commuting) fusion algebra 



E N llt2 h n l3 . (2.22) 



By fl2.18|) , rii = I, the unit matrix, and by (|2.17|) 



rij* = nj '. 



Conversely, given a set of matrices with non-negative integer elements, satisfying rii* = 
nj, ri\ = I and the fusion algebra, they form a commuting set, and thus each rii commutes 
with its transpose. These matrices are thus normal matrices that may be diagonalized in an 
orthonormal basis. Their eigenvalues are of the form S^/ 'Sij for some j, and they may thus 
be written in the form (|2.16|) . If one pretends to determine the spectrum £ from the n's, one 



has to impose also that j = 1 appears only once in £, as a manifestation of the uniqueness 
of the vacuum. 

We thus conclude that the search for orthonormal and complete solutions to the Cardy 
equation is equivalent to the search for N valued representations of the fusion algebra satis- 
fying nj = rii*. 

This is the first important result of this paper, already presented succinctly in ||32|| . The 
fact that some solutions to the Cardy equation were associated with representations of the 
fusion algebra had been noticed before. In his seminal paper |H| , Cardy considered the case 
of "diagonal theories" (for which £ = X) and showed that the matrices were nothing other 
than the fusion matrices Ni, thus obtaining an alternative and more intuitive derivation of 
the Verlinde formula. In an antecedent work by Saleur and Bauer ||, other solutions had 
been obtained in non-diagonal theories, starting from their lattice realization, and the fact 
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that these rii coefficients satisfied the fusion rules had been emphasized in [BUf . More recently 



Pradisi, Sagnotti and Stanev |14|, [15], |30[ proposed a different argument to the same effect, 
where a notion of completeness of boundary conditions is also playing a crucial role. 

Solutions such that all matrices rij may be written as rii = (ni)j © (72.2)1 after the 
same suitable permutation of rows and columns can be called reducible. They describe 
sets of decoupled boundary conditions. We thus restrict our attention to irreducible sets of 
matrices. 



2.3 WZW sl(2) theories 

For theories with the affine (current) algebra sl(2) as a chiral algebra, the problem of classi- 
fying representations of the fusion algebra was solved long ago |39| . The integrable highest 
weight representations of s/(2)fc at level k G N are labelled by an integer 1 < j ' < k + 1, 
Sij = J sin and the Cardy equation says that the generator n 2 = n 2 * has eigenvalues 



k+2 



q 

= 2 cos The only symmetric irreducible matrices with non-negative integer entries 



and eigenvalues less than 2 are the adjacency matrices of A-D-E-T graphs |^0[ of Figure 2 
(see also Table 1). The "tadpole" graphs are given by T n := A2 n /%2- Here the level k is 
related to the Coxeter number by g = k + 2. Only the A-D-E solutions are retained as their 
spectrum matches the spectrum of sl(2) theories, known by their modular invariant partition 



functions [41] . For a theory classified by a Dynkin diagram G of A-D-E type, the set £ is 
the set of Coxeter exponents of G as in Table 1. The matrices n« are then defined recur- 
sively by rii = I , n 2 = G and by equation ( |2.22| ) which reduces here to n i+ i = n 2 rii — rij_i, 
i = 2,3, ... ,k. They are the well known "fused adjacency matrices" or "intertwiners" Vi, 
studied in and whose properties are recalled in Appendix B. One verifies that all 

their entries are non-negative integers. This set of complete orthonormal solutions of the 
Cardy equation for s/(2) theories is unique up to a relabelling of the states \a). 



2.4 Minimal s/(2) models 

The classification of c < 1 minimal models in the bulk is given in terms of a pair of Dynkin 



diagrams (A, G) where G is of A-D-E type [|1]. Let h be the Coxeter number of Ah-\ and 
g the Coxeter number of G as given in Table 1. Then the complete A-D-E classification is 

( M(A h _i, A g -x) 



M(A, G) 



M(A h _i, D 



(S+2)/2y 



g even 



M(A h ^,E 6 ) 
M{A h ^,E 7 ) 
VM{A h „ u Ez) 

with h, g > 2 and central charges given by 

Q(h-gf 



(2.23) 



1 - 



hg 



(2.24) 
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Figure 2: The A-D-E-T graphs. 
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m e Exp(G) 


A n 


n + 1 


1,2,3, ...,n 


D n 


2n-2 


1, 3, 5, . . . , 2n — 3, n — 1 


Eq 


12 


1,4,5,7,8,11 


E 7 


18 


1,5,7,9,11,13,17 


E 8 


30 


1,7,11,13,17, 19,23,29 


T 


2n + 1 


l,3,5,...,2n- 1 



Table 1: The Coxeter number g and Coxeter exponents m of the A-D-E-T graphs. 
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We will use M. (A, G) to denote these minimal theories. Since g and h must be coprime and 
g is even for all non-A cases, one may always assume, at the price of a possible interchange 
in the (A, A) case, that h is odd, h = 2p + 1. 

Some members of these series are identified as follows: 

M(A 2 , A 3 ) = critical Ising c = 1/2 

M(A 4 , A 3 ) = tricritical Ising c = 7/10 

M(Ai, D 4 ) = critical 3-state Potts c = 4/5 (2.25) 

A4(Aq, D4) = tricritical 3-state Potts c = 6/7 

We will use G to denote both the Dynkin diagram and its adjacency matrix. We use 
r, r 1; r 2 to denote nodes or exponents of A h _i, s, Si, s 2 for the nodes (or exponents) of A g _i, 
b for the nodes of G. We refer the reader to Appendix B for more data on these 
matrices and their eigenvectors. 

If Exp(G) denotes the set of exponents of G (see Table 1), the modular invariant partition 
function of Ai(A h _i, G) reads 

j h-l 

Z = -^ \Xrs(q)\ 2 + off-diagonal terms . (2.26) 

r=l seExp(G) 

The factor | removes the double counting due to the well-known identification of the (r, s) 
and (h — r, g — s) representations of the Virasoro algebra. The diagonal terms in Z, i.e. the 
left-right symmetric (highest weight) states in the spectrum are thus labelled by the set 

£ = {j = ( r , s) = (h - r, g - s); 1 < r < h - 1; s e Exp(G)} . (2.27) 

Each of the unitary minimal models A4(Ah_i,G) with g — h = ±1 can be realized 
as the continuum scaling limit of an integrable two-dimensional lattice model at criticality, 
with heights living on the nodes of the graph G. In particular, the critical series with 
g — h — 1 is associated with the A-D-E lattice models |^3l [0J and the tricritical series with 



g — h = — 1 is associated with the dilute lattice models [45, [46|]. In the non-unitary cases 



the associated lattice models [[47| possess negative Boltzmann weights. In the construction 
of the corresponding lattice models as well as in the description of boundary conditions, it 
turns out that the two diagrams of the pair (A, G) do not play a symmetric role. 

According to the discussion of the previous section, we have to study the fusion algebras 
of minimal models and their (integer-valued) representations. The Verlinde formula for the 
fusion coefficients makes use of the matrix S of modular transformations of characters 

S rsr , s , = /A(_i)(^)(^') sin7rrr 'i^ sinnss '^A (2.28) 
V 9h h g 

with the restriction r,r' odd (or any equivalent condition). The fusion coefficients are then 
found to be tensor products of those relative to the si (2) algebras of level h — 2 and g — 2, 
up to a symmetrization which accounts for the identification (r, s) = (h — r, g — s) 

Ks,r>s' r " s " = N T /'N S /' + N rr , h - r "N ss , 9 - s " . (2.29) 
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This may be regarded as the regular representation of the fusion matrices Af rs of the Virasoro 
algebra of central charge ( 2.24Q . Our problem is to find the general non-negative integer 
valued representations of this algebra. One observes that Af rs = Ari-A/L and that the algebra 
is thus generated by A/21 and Afu- Also, the eigenvalues of A/12 and A/21 are of the form 

Sl2 ' r ' s ' (-1)^2 00s tts'^ - = (-IT' 2 costs'- , (2.30) 



ll,rV 



S 

^21,r's' 1 -i\s' rt 1 9 /o Ol\ 

— — — = (—1) 2cos7rr— . (2-31) 
e>n,r's' h 

with 1 < r' < h — 1, 1 < s' < # — 1 and again (r', s') = (h — r',g — s'). 

Turning now to a general (integer valued) representation n rs of the fusion algebra, it 
is still true that it is generated by 7112 and 7121. In addition, we want the spectrum of the 
n rs to be specified by the set of "exponents" £ of ( |2.23| ), that is (r', s') in ( |2.30| - |2.31] ), with 
the eigenvalues labelled by s' appearing with some multiplicity in general. To remove the 
redundancy in the labelling of eigenvalues, we will usually take r' odd, r' = 1,3, ■ ■ ■ , h — 2, 
and (s', a) G Exp(G). In the sequel, we will drop this explicit notation for multiplicities. We 
know of course a solution to this problem, namely 

n TS = N r ®V s + N h ^ r Vg- S (2.32) 

in terms of the fusion matrices N of si (2) at level h — 2 and of the intertwiners V of type G 
introduced in the previous subsection (see also Appendix B). More explicitly, this describes 
a solution to the Cardy equation between boundary states (n, a) and (r 2 , b) 



n 



rs;(r\ ,a) 



(r2,6) = N rri r Wj + N h _ rri r2 V g _ sa b , (2.33) 



with 1 < r, ri, T2 < h — 1 = 2p, 1 < s < g — 1, and a, b running over the nodes of the Dynkin 
diagram G. 

Because of the properties of the N and V matrices recalled in Appendix B, it is readily 
seen that 

„ (r 2 ,b) (h-r 2 ,j(b)) _ „ ,,(^2,6) (n VA) 

for an automorphism 7 acting on the nodes of the graph G: this is the identity except for the 
A -Dodd and E'g cases, for which it is the natural Z 2 symmetry of the diagram. We conclude 
that this solution describes boundary states of M.(A h _i, G) labelled by pairs (r, a) of nodes 
of the Ah-\ and of the G graph, with the identification 

(r,a) = (/i-r j7 (a)). (2.35) 
One checks that the number of independent boundary states |(r, a)) is 

number of independent boundary states = —(h — l)n (2.36) 
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with n the number of nodes of G, or the number of its exponents. This number (2~. 



coincides with the number of independent left-right symmetric highest weight states |r, s) ® 
|r, s) in the spectrum of the theory on a cylinder, i.e. with the cardinality of the set S, as it 
should. 

With such boundary states, the cylinder partition function reads 



'(ri,a)|(r 2 ,b) 



J (r 2 ,b)\( ri ,a) 



(2.37) 



Here the sum runs over 1 < r < h — 1, I < s < g — 1. 

Let us look more closely at (|2.33|) . There exists a basis in which (|2.33| ) takes a factorized 
form. Indeed one may use the identifications (r, s) = (h — r, g — s) and ( [2.35D to restrict 
r, ri,r 2 to odd values (recall that h = 2p + l is odd). Then Nh- rri r2 = 0, the r.h.s. of ( |2.33| ) 
factorizes and the following expressions 



(r"s") 
(r,a) 



V2S rr »ipi', r, r" odd, s" G Exp(G) 



(2.38) 



written in terms of the modular matrix of si (2) at level h — 2 and of eigenvectors ■?/> of G, 
are readily seen to be eigenvectors of n rs . Their eigenvalue is of the form S rs y s > /Sny s / after 
some reshuffling r", s" — * r',s'. 

One also shows (see Appendix C) that there exists a basis in which 



ni2 



I P ®G 



n 2 i 



V 

r 

V 



G 



r 

o 



r 



g; 
\ 

r 



(2.39) 



(2.40) 



in terms of the tadpole T p adjacency matrix and of T, the matrix that realizes the automor- 
phism 7: T a b = S a7( b) • 

Conversely, suppose we only know that the representation n rs has a spectrum specified 
by the set of exponents £. The question is: are these spectral data sufficient to guarantee 
that the only n rs are of the form ( |2.33j ) in a certain basis? A proof of this fact is relegated to 
Appendix C. Notice that our discussion has assumed the classification of modular invariants 
to be known. It should be possible to extend it as in the case of WZW si (2) models and to 
classify the representations of the Verlinde algebra without this information. A few spurious 
cases involving tadpoles etc. would then have to be discarded. 

To recapitulate, we have proved that the only representations of the fusion algebra of 
minimal models are given by ( |2.33| ). To our knowledge, this is the first proof of the uniqueness 
of these (complete orthonormal) boundary states of minimal models. 
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Some physical intuition about the meaning of these boundary conditions may be helpful. 
For this we appeal to the lattice realization of the minimal model as a generalized height 



model on the graph G (see flSlfl ). A boundary condition of the type (1, a) describes a fixed 
boundary condition, where the height of the model is fixed to value a on the graph. The 
interpretation of the other label r is less intuitive. The boundary condition (r, a) is realized 
by attaching an r-times fused weight to height a. 

The expression ( |2.37| ) for the cylinder partition function encompasses and generalizes 
cases that were already known: 

• From the work of Saleur and Bauer |5] who discussed boundary conditions in lattice 
height models of A-D-E type on a cylinder in which the heights on the boundaries 
are fixed to the values a respectively b. They showed that in the continuum limit, the 
partition function reads 

Zb\a = ^2 VsaXls ■ 
s 

• From the work of Cardy |TD[ who showed how to construct new boundary conditions 
by fusion. 



From the work of Pasquier and Saleur [p9|| , who interpreted the pair of relations 



Z \tm\'S ] = ExuVj (2.41) 
^(rfr = Xu (2-42) 

as expressing the decomposition of the representation of the Temperley-Lieb algebra 
on the space of paths from a to b on graph G onto the irreducible ones on the paths 
from 1 to s on graph A g _i, see point (ii) at the end of Appendix B. 

Examples 

Let us illustrate these expressions of boundary states by a few simple cases. In the Ising 
model (the (A 2 , A 3 ) minimal model), h = 3, G = A 3 , thus n = 3 and there are |(3 — 1) x 3 = 3 
boundary states, generally denoted [[Kj] +,— and /. On the lattice, the first two describe 
fixed boundary conditions on the spin a = 1 or — 1 respectively, while / corresponds to free 
boundary conditions. 



It is then instructive to consider two related examples, see also |22|, |28|. The first is 
the c = 2 solution of sZ(2) 4 at level 4, and the other is its cousin, the c = 4/5 minimal 
(3-state Potts) model, already mentioned in Section 2.4 and labelled by the pair (At,-D 4 ). 
In the former case, we find four boundary states, labelled by 1 to 4, that we attach to 
the nodes of the -D4 diagram. All these states satisfy the required boundary conditions. 
The set of exponents is £ = {1,3,3,5}. But this -D4 s/(2) model is also known to result 
from the conformal embedding of s/(2) 4 into s/(3)i. Regarded as an si (3) theory, the model 
admits three boundary states satisfying the more restrictive sl{3) conditions (L n — L_ n )\a) = 
(J n + f2J_ n )|a) = 0, where the choice of fl corresponds to the diagonal set (J,j). These 
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three boundary states may be regarded as the three nodes of a triangular graph A^ (see 
Appendix D and Figure 10), or as the three extremal nodes of the -D4 diagram that have 
survived the additional si (3) constraint. 

The discussion of the Potts model is quite parallel. From the minimal model standpoint, 



it is the (A 4 , -D4) model, h = 5, n = 4 and there are 8 boundary states [22. 28 



A = (1, 1) = (4, 1), B = (1, 3) = (4, 3), C= (1, 4) = (4, 4) 
BC=(2,1) = (3,1), AC = (2, 3) = (3, 3), AB = (2,4) = (3,4) (2.43) 
ABC =(1,2) = (4, 2), iV=(2,2) = (3,2) 

On the lattice, the first three A, B, C describe fixed boundary conditions where the "spin" 
takes at each site of the boundary one of the three possible values. The mixed boundary 
conditions AB, BC, AC describe boundary conditions where the spin on the boundary can 
take on two values independently. The boundary conditions ABC and iV are free boundary 
conditions but for N the weights depend on whether adjacent spins are equal or not. 

The model may also be regarded as the simplest W3 model. In that picture, one may 
impose more stringent boundary conditions. Only the six states denoted above A, B, C, 
AB, BC, AC satisfy the additional condition (W^ ] +fW( 3 )_„)|a) = 0. They correspond to 
the extremal nodes of the pair (T 2 ,^4) or, alternatively, to the nodes of the pair (T 2 , A^ 4 -*). 
As will be discussed in more detail in Section 3, the subset of these nodes, to be denoted 
T, can be identified in both examples with the representation labels of the corresponding 
extended chiral algebra. The matrix elements ijj J a for a G T satisfy |53j the relation 



s 



cxt 

°di2 



Qext 



where {j} denotes the orbit of the exponent j with respect to the Z 2 automorphism 



and Suffix is the modular matrix of the extended theory. This relation implies that \a) = 

gext 

^0'} / — i- e -5 we can identify J2je{j} = 10})) w ith an extended Ishibashi 

state. The missing boundary condition corresponds to a twisted boundary condition from 
the point of view of the extended algebra. 

We conclude that, as expected, the number and nature of the boundary states reflect 
the precise conditions that they are supposed to satisfy. 



3 . Graph Fusion Algebras 

According to the discussion of Section 2, given a certain chiral algebra A, the sets of complete 
orthonormal boundary states of RCFTs consistent with this algebra are classified by repre- 
sentations of the Verlinde algebra of A, on matrices rij with non-negative integral entries. 
Or stated differently: given a certain RCFT with a chiral algebra A, the sets of complete 
orthonormal boundary states of this theory consistent with this algebra are classified by rep- 
resentations of the Verlinde algebra of A, on matrices with non-negative integral entries, 
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with eigenvalues specified by the diagonal part 8 of the spectrum. These matrices may thus 
be regarded as the adjacency matrices of a collection of \I\ graphs. In practice, it is sufficient 
to look at the smaller number of matrices representing the generators of the fusion ring. For 
example, one matrix in the case of sl(2) considered above, or the N — 1 matrices associated 
with the fundamental representations, in the case of sl{N). 

The simplest case is given by the regular representation of the fusion algebra, when 
the matrices rii are the Verlinde matrices themselves, rij = Aj. This is the case of so- 
called diagonal theories, when all representations of the set X appear once in the spectrum 
Spec = G X}. This may be regarded as the case of reference from several points 

of view: it was the first case analysed in detail ||10|| ; the corresponding graphs are playing 
a central role; and finally in that case, Cardy was able to provide a physical argument 
explaining why the fusion matrices arise naturally. It is the purpose of this section to extend 
these considerations to more general solutions. We shall find that the role of the fusion 
matrices in the arguments of Cardy is now played by two sets of matrices. The first is the 
set of matrices that describe the coefficients of the cylinder partition function; the second 
is a new set of matrices N a , forming what is called the graph fusion algebra. 

On the other hand, since we know that the cylinder partition functions, or equivalently 
the matrices rij, contain some information about the bulk theory, through the knowledge 
of the diagonal spectrum £, it is expected that this classification of boundary conditions 
should have some bearing on the classification of bulk theories, namely on the classification 
of torus partition functions and on bulk structure constants. Remarkably, this programme 
works even better than expected and the two classification problems seem to be essentially 
equivalent, at least for type I theories (see end of Section 1). This will be explained in 
Section 3.3 below. 



3.1 More on graphs and intertwiners 

Suppose we have found a solution to the Cardy equation, namely a set of n x n matrices 
( |2.16| ), (rii) a b , i — 1, • • • , \X\, a, b — 1, • • • , n. What was said in detail in Section 2 and in 
Appendix B in the case of sl(2) can be repeated here. As their entries are non-negative 
integers, these matrices may be regarded as adjacency matrices of a set of |X| graphs G i: 
with n = \Gi\ = \G\ nodes. We shall refer collectively to these \I\ graphs as "the graph 
G" , whereas the basic solution provided by the A's themselves will be called "the A graph" , 
(borrowing the notation from the s/(2) case). The eigenvalues of the matrices nj are specified 
by a set Exp(G) in the sense that they are of the form Sij/Sy, (j,a) G Exp(G). Moreover 
Exp(G) = £ if the RCFT is given and the diagonal spectrum 8 is known. But in general, 
the determination of the set Exp(G) is part of the problem. The fundamental relation 
T,b( n i)a b {nj)b c = Efc Nij k (n k ) a c may be interpreted in two ways: 

• Regarded as \G\ x \G\ matrices, the matrices rii form a representation of the fusion 
algebra ( g22| ). 
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• Regarded as a \X\ x \G\ rectangular matrix, each matrix n a for a fixed, {n a )j b '■= rij a h 
intertwines the representatives iV, and rij in the two representations A^n a = n a nj, or 
more explicitly 

E%V = E%V. (3.1) 

k b 

We shall thus occasionally refer to the matrices n.j as "intertwiners" . 

The case of graphs and intertwiners pertaining to si (2) theories has been discussed at 
length in Section 2 and in Appendix B. In Appendix D, we outline the discussion of si (3). 
In that case, the fusion algebra is generated by two matrices ri(2,i) and ri(i,2) (labelled by the 
two fundamental (shifted) weights of sl(3)), and as these two representations are complex 
conjugate to one another, the matrices are related by transposition n^i) — nT 12 y Then 



according to ( |2.22|) , is given by the same polynomial of ^(1,2) and ^(2,1) with integral 
coefficients as that representing Ni in terms of Nn^) and Npu. It is thus sufficient to list 
all possible graphs representing the matrix Ji(2,i), provided one checks that all rii have non- 
negative integral entries. In contrast with the case of si (2), no complete solution is known 
for sl(3). The current state of the art is presented in Appendix D with tables, figures and 
relevant comments. 

3.2 Graph fusion algebras 

To see what is playing the role of the fusion algebra in the argument of Cardy, we have 
to introduce the graph fusion algebra. The graph fusion algebra, as first discussed by 
Pasquier IpTfl , is a fusion-like algebra attached to a connected graph G. Let i/; 3 a be the 



common orthonormal eigenvectors of the adjacency matrices G labelled by j G Exp(G). In 
general, these eigenvectors can be complex. In the case of degenerate eigenvalues the associ- 
ated eigenvectors need to be suitably chosen. We assume that the graph has a distinguished 
node labelled 1 = 1* such that > 0, for all j e Exp(G). 
One then defines the numbers 

E mm: (3. 2) 

jeExp(G) W\ 

and the matrices N a with elements (N a ) b c = N ab c satisfy N ab c = N ac * b * and iVj = A^ a *. 
Because of ort honor mality, N± = I. Since each matrix N a has a single non- vanishing entry in 
the row labelled 1, namely (N a )i b = (Ni) a b = 5 a b, the matrices N a are linearly independent. 
The N ao c are the structure constants of the graph fusion algebra satisfied by the iV matrices 

N a N b = Y,N ab c N c (3.3) 

c 

which is an associative and commutative algebra. Of course, if the graph G is of type A, 
this boils down to the ordinary Verlinde fusion algebra since the matrix ip of eigenvectors is 
nothing but the modular matrix S. 
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Since the N a and rii matrices have the same eigenvectors, it is easy to derive the matrix 
relation 

*N a = J2n ia b N b . (3.4) 



nil 

b 



In particular for a — 1, Ni — I, and all n,i appear as linear combinations with non-negative 
integer coefficients of the iV's 

nJ = ^2^i c N ca b (3.5) 



Alternatively ( |3.4|) may be used as a starting point to reconstruct the graph algebra, as 
explained in Section 3.5. 

It should be stressed that the definition of a graph fusion algebra is not unique. In 
general, it depends on the choice of the distinguished node 1 and, when there are degenerate 
eigenvalues, also on the choice of the eigenvectors ip J a . To view the graph fusion algebra as a 
proper fusion algebra we would like the structure constants N a b c to be non-negative integers. 
But even the rationality of these numbers is not obvious and it is therefore surprising that, 
for appropriate choices of the f/>'s and of node 1 and for most cases, they turn out to be 
integers of either sign. Among all the examples known to us in sl(N) theories, 2 < iV < 5, 

(12) 

it fails in only two cases: the graph called in Figure 11, for which there is no node 1 
satisfying 1 = 1*, and whose N algebra involves fractions of denominator 4; and a graph in 
the s/(4) 4 theory, f53fl , in which half-integer N ab c of either sign occur. Adopting ( |2.22|) , 



in the framework of subfactors theory the latter example has been reinterpreted by Xu 
by trading commutativity of the graph fusion algebra for integrality. 

Finally the non-negativity of the N is only possible for certain graphs which we call 
proper fusion graphs. For example, for the sl(2) theories, the A-D-E graphs that admit a 
proper graph fusion algebra are 

proper A-D-E graphs = A n , D 2q , E 6 , E 8 . (3.6) 

The choice of distinguished node for the sl(2) A-D-E graphs is explained in Appendix B. 

We note that the set of proper sl{2) fusion graphs matches the modular invariant parti- 
tion functions listed as "type I" at the end of Section 1. The situation is somewhat different 
for sl(3) graphs, for which we have to introduce a further distinction. In this case, some 
graphs with non-negative iV's are not associated with type I theories. We reserve the termi- 
nology "type I graph" for those graphs associated with type I theories (see Appendix D and 
Tables). Moreover, as is clear from the 5/(4)4 example above, some type I modular invariant 
partition functions are associated with graphs with non- integer and/or non-positive N a b c - In 
the following, we discard these exceptional cases and restrict ourselves to type I graphs that 
are associated with type I RCFTs. In general, the question of precisely which graphs admit 
type I fusion algebras should be related to the classification of type I RCFTs, and thus is a 
very interesting open question. 
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3.3 Fusion rules and block characters 



Given a solution to the Cardy equation, that is a set of partition functions 

z a\b(q) = J2 n ™ b Xi((l) (3.7) 

and the corresponding graphs, we assume as in Section 3.2 above that there exists a special 
node called 1. We then introduce the combinations of characters (or "block characters") 

Xc(q) = ^hicxM) ( 3 - 8 ) 

where 

n ic = nn c (3.9) 
is referred to as the basic intertwiner. Thanks to (|3.5| ), ( |3.7|) may be rewritten as 

Z a \ b {q) = Y,tica b xM) (3-10) 

c 

where the coefficients are now given by the structure constants of the graph fusion algebra 
of G. 

Equation ( p.l(J| ) is a mathematical identity and as such is valid and consistent indepen- 
dent of the choice of the distinguished node and eigenvectors of G. Physically, however, the 
case where the N a b c are non-negative integers is the most interesting. In that case, following 
Cardy's discussion [jnj, it is suggested that N ab c gives the number of times that the propa- 



gating mode or representation c appears in the strip or cylinder with boundary conditions 
a and b. Thus if G is a type I graph, i.e. if the structure constants N ab c are non-negative 
integers, we have a possible interpretation: the nodes a of the graph(s) G under considera- 
tion label a class of representations of some extended chiral algebra. The blocks Xa are their 
characters, and the integer coefficients N a b c are their fusion coefficients. 

To probe this interpretation, let us see how it confronts the results "in the bulk", in 
particular how it is consistent with the form of the torus partition function. There, it has 
been observed already long ago that (for type I theories) the torus partition function (cf 
( |1.5| )) may be recast in the form 



^torus — X/ \Xa\ 2 i (3-11) 
aeT 



i.e. as a diagonal sum over a subset T of block characters. The subset T corresponds 
to a subalgebra of the N algebra, in the sense that if a, b G T, N a b c ^ only if c £ T. 
This interpretation of multiplicity of representation i in the block c, that was first 

observed empirically |5l| , was subsequently derived in a variety of cases of type I sl(N) 



theories based either on conformal embeddings or on orbifolding p3| , |54fl . More recently, it 
appeared as an important ingredient in the investigation of the algebraic structure underlying 
these theories 0, |5T|j . 
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The following interpretation is thus suggested. The nodes a G T label representations 
of the maximally extended algebra A' of the RCFT (of type I) under consideration. The 
subalgebra of the N algebra is the conventional fusion algebra of this RCFT. The other nodes 
a ^ T might label other "twisted" representations. The entire N algebra would describe the 
fusion of all, twisted and untwisted, representations of A'. This interpretation in terms of 
twisted representations seems corroborated by the fact that some \ a are known also to occur 
in partition functions on a torus in the presence of twisted boundary conditions. The fact 
that general boundary conditions on a cylinder also appeal to these representations was first 
observed in the Potts model in p2[ . See for some work in this direction [jSBfl , and the more 
systematic developments |29| along the lines of f63| . The concept of twisted representations 
in other cases, like conformal embeddings, remains to be understood. 



Having discussed the situation for type I theories and graphs, we return to RCFTs of type 
II. There the situation is more elusive. On the one hand, as discussed above, the boundary 
conditions on a cylinder are labelled by nodes of an improper graph G, and although we can 
still write an expression of the form ( p.lOj ), its physical interpretation is unclear. On the 
other hand, from ( |1.6|) , we know [|73| that the torus partition function may be expressed in 
terms of block characters pertaining to a "parent" type I theory with graph G' 

Ztoms = E Xa(q) (Xc(a)(q)) , (3-12) 
aeT 

where T is once again a subset of the nodes of G' corresponding to a subalgebra of the N 
algebra, and ( is an automorphism of that subalgebra -Nf(a)c(6)^ = N a b c - 

For example, the s/(2)i6 theory labelled by the Dynkin diagram Ej is known to be 
related in that way to the D w theory. Their respective torus partition function read 

IXi + Xn\ 2 + |X3 + Xisl 2 + |X5 + X13I 2 + |X7 + Xnl 2 + 2 1 X9 1 2 (3-13) 

E lE4 Dl0) X,| 2 (3.14) 

0=1, 3,5, 7,9,10 i 

aeD 1Q 

IXi + X17I 2 + |X5 + X13I 2 + |X7 + Xnl 2 + IX9I 2 + ((X3 + Xis)Xg + c - c -) (3-15) 
E (EC 10 ^) (E^xX , (3.16) 

0=1,3,5,7,9,10 \ j / \ i / 

aeD 10 

with ( exchanging the two nodes 3 and 10 of the D\q diagram. 

It seems that the parent graph G' also plays a role for cylinder partition functions of 
type II theories. Indeed, to obtain cylinder partition functions expanded with non-negative 
coefficients in terms of block characters, we just have to expand in the block characters of 
G' . Specifically, we find 

= E n£ G ' )b xf Hi) 

c&G' 



7 {Dio) 
•^torus 



7 (E T ) 

•^torus 
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where the G-G' intertwiners are given by 

„(GG')b _ V- 1pf )m ,AG)m ,{G) m * 
mGExp(G)^l m 

These turn out to be non-negative integers n^ G '^ b > and satisfy the G' graph algebra. 
Here it is assumed that the distinct exponents of Exp(G) are in Exp(G') and that the sum is 
over exponents of G counting multiplicities. Moreover, if there is more than one eigenvector 
of G' corresponding torn £ Exp(G), then any of these eigenvectors can be matched with the 
given m G Exp(G). The formula can be derived in the same way as our previous formulas. In 
particular, this formula applies for G = E-j and G' = D w . In this case there is an ambiguity 
as to which D\q eigenvector is taken for m = 9 but in fact one can take either. The matrices 
are changed by the Z-i symmetry but the cylinder partition functions agree. The formula 
also holds for the type II sl(3) theories G = £ 2 (12) , £ 4 (12) and £ 5 (12) and G' = V^ 12 \ 




L L 

Figure 3: The torus partition function reconstructed from two cylinder partition functions. 

Putting everything together, we finally observe that in general for a "rectangular" torus 
with two periods 2L and iT, made by pasting together two cylinders, (see Figure |3|), 

7(G) _ y(G') y{G') /q 1 7 \ 

Aotus - 2^ Z a|l Z C(a)|l ' \ 6 - L <) 
aeT 

i.e. the partition function may be obtained as the sum over a special set of boundary 
conditions of cylinder partition functions. This expression is of course deeply rooted in all 
the connections between bulk and boundary theories, open and closed strings, etc, but still 
we find its simplicity intriguing. 



3.4 Examples 

More examples can be given to the previous general scheme. 

• sl(N) : the classification of the representations of the fusion algebra of sZ(iV)fc is a 
well posed problem on which we have only partial results. In particular, classes of 
graphs pertaining to si (3) as well as some cases for higher N have been expounded 



from various standpoints in |39|, |51], |53], |54], p9|] , (see Appendix D). In all known 
cases, the previous discussion may be repeated: intertwiners, type I graphs, and other 
concepts introduced above, still apply. We refer the reader to the above references. 
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The case of sl(N)i may be described in detail. The representations of sl(N)i are 
labelled by an integer < % < N — 1 (we depart here from our previous convention, 
with i = denoting the identity). The fusion rules are isomorphic to the addition of 
integers modulo N, and the algebra defined by NiNj = N i+ j moi iN is thus generated by 
N\, Nj = (Ni)i . The eigenvalues of Nj are exp2injl/N. The regular representation 
is provided by N x N matrices, generated by (Ni)^ = 5j,i+i m odiv- All the previous 
eigenvalues are reached once and we may thus say that the system has "exponents" 
I = 0, • • • , N — 1 (mod N). In general, a representation {rii} of the fusion algebra is 
associated with each divisor q of N = p.q, including q = 1 and q = N: q denotes the 
order of the matrix ri\ which is qxq dimensional and such that (jii) a b = Sb- a ,imodq, for a 
labelling of the nodes a,b = 1, • • • , q. (This exhausts all integer-valued representations 
of the algebra. Indeed the conditions that nj = njv_i and nitiM-i = I imply that 
the only entries of n\ are and 1, and that ri\ is a permutation matrix. Being of 
order q and indecomposable, n\ is a matrix of a cycle of length q. Q.E.D.) Obviously 
the matrices rii = (ni) ? are all integer-valued, n = I, and nf = rii* = n^_i. The 
graph of adjacency matrix n\ is an oriented g-gon. In that case, we may say that the 
q exponents are 0,p, ■ • ■ ,p(q — 1). 

This census of representations of the sl(N)i fusion algebra matches almost perfectly 
that of modular invariant partition functions carried out by Itzykson |6(| and Degio- 
|6l . We recall that according to these authors, a different modular invariant is 
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associated with each divisor of n, where n = N if N is odd and n = N/2 if it is even. 
Thus, only the case N even, q — 1 has to be discarded in our list of representations of 
the fusion algebra, as it does not correspond to a modular invariant. 



3.5 More on graph algebras 

In Section 3.2, we have introduced the matrices iV by (|3.2| ) and derived ( |3.4|) . Instead of 
looking at the graph as a collection of points we can look at it as a collection of matrices N, 
providing a basis of a commutative, associative algebra with identity, and an action of the 
intertwiners rii given by (|3.4|) , that is, we take ( |3.4| ) as a starting point. Given the graph G, 
in particular the coefficients rin c , it is possible in many cases to invert ( |3.5| ) and solve for 
N a as linear combinations with integral coefficients of the intertwiners rii, or equivalently, as 
polynomials of the fundamental adjacency matrices. Similarly the relation fl3.4p written in 
terms of the eigenvalues 7 3 -(i) = fj 1 , = 

7j(«) 7o(i) = Yl n i* W) ' * G 8 ' j G J ' ( 3 - 18 ) 

b 

is a recursive relation determining (the rows of) the eigenvector matrix ij} % a . In general, 
typically in the presence of degenerate eigenvalues, the matrix ip l a is not determined uniquely, 
or alternatively, ( |3.5| ) cannot be inverted for all a G V. For the type I cases, however, as 
explained above, there exists an extended fusion algebra isomorphic to a subalgebra of the 
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graph algebra, so that the extended fusion matrices can be identified with a subset 



iV a with the nodes a E T C V ||51|| . In most of these type I cases one can solve for all N a 
in terms of the n^s and N^ 1 , or alternatively express ip J a in terms of the modular matrices 
Sij , S^ Bi . A particularly simple subclass of Type I for which one can go quite far in the 
programme of reconstructing the graph G and all the related structures is presented by the 
orbifold theories, in particular the ones associated with groups generated by simple currents. 
In our approach they can be described by graphs obtained by 'orbifolding' the fundamental 
graphs of the initial diagonal theory, the simplest example being provided by the WZW sl{2) 
D 2 i series obtained by "orbifolding" the Dynkin diagram over the Z2 group generated 

by the automorphism 7. In these cases as well as in their sl(N) generalisations defined in 
|62"H involving the group Z^v, one can algorithmically construct the eigenvector matrix, see 
Appendix B for an illustration in the simplest N = 2 case. In a different approach, using 
tools similar to the original orbifold treatment of |S3], an elegant general formula for the 



eigenvector matrix was derived recently in p9| . It should be noted that the same graph 



(orbi)folding procedure leads also to type II graphs, e.g., the sl(2) D odd series, or their s/(3) 
generalisations for k 7^ mod 3, see Appendix D. 

We have assumed up to now in this discussion that the graphs are already known. On 
the other hand the relations ( |3.4| ), ( |2.22| ) can be taken as the starting point for finding new 
graphs, typically "exceptional" graphs not covered by the previous orbifold constructions. 
Since any graph in the vicinity of the identity resembles the original "diagonal" (A) graph, 
one can first try to identify n^s for which the r.h.s. of ( |3.5|) reduces to one term, i.e., 
nn a = S aai and hence one can identify n,i = N ai . According to ( |3.18|) this also determines 



tPa- ty once tp{ is known. This is a problem which is reduced to the computation of 
some Verlinde fusion multiplicities. Indeed let us take the first matrix element a = b = 1 of 
the matrix relation Q2.22 ) we have 



E V c = E V = E N u j ( 3 - 19 ) 

C I l£p 

where in the last sum p = {I G T\nn l 7^ 0} and / is counted rin 1 times. Let us assume 
first that in ( 3.19|) i = j. Whenever the r.h.s. of ( |3.19|) is equal to 1, since by definition 



rin a are integers, the l.h.s. summation reduces to one term, i.e., we recover nn a = 5 aai . 
Furthermore plugging this into the l.h.s. of ( |3.19[ ) taken for j 7^ i we recover nji ai as being 
given by the sum of Verlinde fusion multiplicities in the r.h.s. of ( |3.i9j ), i.e., we determine 
the multiplicity with which iV a . appears in rij, see ( |3.5| ). Similarly, a value 2 or 3 for the 
r.h.s. of ( |3.19|) with i = j would lead to 2, respectively 3 terms in (|3.5|), while 4 could be 
interpreted either as leading to 4 terms with multiplicity one, or 1 term with multiplicity 
two, i.e., ni\ a = 25 aai . What we only need in order to check all these possibilities is to 
know the content of the set p, i.e., Tin 1 . This data is provided in type I theories for all of 
which p encodes the content of the identity representation of the extended algebra. More 
generally, nn a = mult £ a (i), identifying a with a representation B a of the extended algebra. 
The relation ( |3.19| ) and its consequences just described are the first steps in a consistent 
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algorithmic procedure proposed by Xu in the abstract framework of subfactors theory 
(see also |56|, 57] for further developments). In particular, the subset of rii which can be 
identified with some N ai are related to "irreducible" sectors with the sum in the l.h.s. of 
( 3.19|) interpreted as a scalar product (ai,aj). The algorithm reduces systematically the 
determination of rii , N a in type I cases to data provided by the Verlinde fusion matrices N 
and N ext . All graphs previously found in [[39], [51], [53|] were recovered in |55j by this method 
and a new example corresponding to the s/(4) 6 modular invariant was found in [Q. 

Finally let us point out that to some extent this algorithm for finding new solutions of 
the equations fl3.4| ) and ( f2.22j ), i.e., new graphs, can be applied to type II cases where we 
do not know a priori Tin 1 , i.e., the set p. One can start with some trial set and compute 
Y, a (w ) 2 — SijJ2i£ P Sij- A first consistency check is that Y^jeei^i) 2 — 1- Then one 
can proceed as in type I. For example the Ej Dynkin diagram may be reconstructed using 
p = {1,9,17} 0. Different (consistent) choices of the set p might lead to the same graph, 
reflecting the possibility of different choices of the identity node. 

In some simple cases it is possible to recover a complete set of boundary conditions by 
applying formula ( |3.4| ) to a known subset in such a way that only one term appears in the 
sum in the r.h.s. In terms of the equivalent formula ( |3.18[ ) for the eigenvalues we obtain a 
new solution ip l b ,b = b(a,j) , by "fusing" a given one ip l a with the Verlinde eigenvalue 7j(i). 
This seems to be the idea of the so called "Affleck fusion conjecture" jl7| , which clearly has 
a restricted application, with the general formulae (|3.4| ), ( |3.18|) being the correct substitute 
for it. 



Another approach to constructing Type I graphs was discussed in [[53] and used to find 
new solutions for higher rank cases. It is based on the use of a relation for the structure 
constants of the Pasquier algebra, the dual of the graph algebra, with structure constants 
labelled by elements in the set Exp(G) and given by a formula analogous to ( |3.2| ), however 
with the summation running over the nodes of the graph; this algebra will be discussed 
further in section 4.4.1 below. 



4 . Bulk and Boundary Operator Algebras 

In this section we investigate the algebras of fields in the presence of boundaries and the 
equations for their structure constants resulting from duality constraints. Our discussion 



parallels that of Cardy and Lewellen ||11|| , but generalises it in two respects: to higher rank 
and to non-diagonal theories. This results in additional multiplicities associated with the 
more general representations of the Verlinde algebra (|2.22|) . Our presentation makes use of 
concepts used by Moore and Seiberg for bulk RCFTs and extends them appropriately for 
this new setting. This leads to a richer structure in the equations and the appearance of 



3 Note added in proof: this was independently discussed in the recent paper |p6[ , see also |73fl . The paper 
[ )S6|| provides a systematic approach in the framework of the subfactors theory to both types of modular 
invariants. 
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a triplet of algebras (n», N a , Mj). Separately these algebras have appeared before but the 
inter-relation between these algebras has not been shown in this context. 

4.1 Ground state degeneracies 



As stressed by Affleck and Ludwig [16], the logarithm of the partition function, in the limit 
L/T —>■ oo, contains not only the universal term proportional to L and to the central charge 
(in unitary theories), but also an L independent term, interpreted as a boundary condition 
dependent "ground state degeneracy" \ng a gb- Indeed in that limit 

In Z b \ a ~ ^ 4tt I + In + - In S u (4.1) 

where as before we denote by 1 the representation of conformal weight 0, corresponding to 
the identity operator. We therefore identify 

9a = ^=- (4-2) 



Thus in unitary minimal models, using ( |2.3S| ) and (|2.28 ), we have the following expres- 
sion for the boundary states 

q(h) is' 

l(r,a)>= E 2i T%r%r |r ' s,)) (43) 

r'odd,s'gExp(G) V 

and their g factor 

_ «l,l|(r,a)) _^ Sjf^l 

9{r,a) ~ , - ^ j ... = , (4.4) 

JqWqW . cW n(g) 

V °n °n V °n °n 

in terms of the modular matrices and of sl(2) at levels h — 2 and <? — 2, |# — /i| = 1. 
For example, for the critical 3-state Potts model, we obtain 

9A = ^ 30^ ^ = °- 550936 (9 a ■ 9ab ■ gABC ■ 9n) = (1 : — : Vs : 1 + ^ Vs) 

in agreement with p2| . As a particular case of ( |4.4j ), the ratio gi r , a )l 9{\,a) equals S±l / , 
in agreement with (A. 3) of ]^2| and the fact that one obtains the boundary state \(r, a)) by 
fusion (in the sense of Cardy) of boundary states (1, a) and (r, 1). 

In non-unitary cases, these expressions have to be slightly amended. If jo denotes the 
representation of smallest conformal weight hj < and assumed to belong to 8, then 

In Z b]a ~ H 4vr ^ + In ^ + In ^ - In S ljo (4.5) 
with c e ff := c — 24hj . Also 

<?a = ^= • (4-6) 
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For simplicity of notation in most of what follows we shall restrict to unitary theories. 
Denoting g a = (l) tt from now on, 



lim Z b \ a e 

L/T->oo 1 



7TC L 

6 T 



7^=(l)a- (4.7) 

One can consider furthermore the partition function with some field insertions at the same 



limit [11 1, [25 1, [26]; we shall normalise them similarly so that only a dependence on g a is 



retained, i.e., (l) n will coincide with the 1-point function of the identity operator in this 
limit. 



4.2 Bulk and boundary fields, OPE 
4.2.1 Boundary fields 



According to Cardy |TD[, boundary conditions can be interpreted as created by the insertion 
of fields b ^j f3 (x) living on the boundary, Im z = 0, x =Re z of the upper half-plane z G H + . 
Here j G X , a,b G V , and (3 accounts for the multiplicity n J0 6 of such fields, to be called 
"of type" Q o J . Thus f3 can be interpreted as a "coupling" index (3 = 1, 2 . . . , rij a b and the 
boundary fields as kind of "chiral vertex operators" (CVO) associated with a second type 
of couplings y a J , a, b G V,j G X. This is a formal analogy since the boundary states 
| a), | b) labelled by a, b are superpositions of Ishibashi states. The multiplicity index (3 is 
traditionally omitted, but it should be stressed that even in the si (2) case, in all but the 
diagonal cases {£ — X), there are always some non-trivial multiplicities rij a b > 1, so most of 
the time we shall retain this index. Since n la 6 = 5 ab the index associated with the coupling 
(J 6 ) takes only one value and will be denoted (3 = lb, or just 1, or, altogether omitted. 
On the other hand in all non-diagonal cases there is a non-trivial subset , j G p} of 
boundary fields, with 1 G p C X, where the set p has been introduced in Section 3.5. 

To make contact with Section 2.2, consider the finite strip w = — log z equipped with 
the Hamiltonian H a b = j (Lq 11 ^ — ^). The space of states is generated by all the descendent 
states created for fixed a,b from the (properly normalised) vacuum state 6 ^ r " j3 (0)|0) by 
the modes of the Virasoro algebra generating the real analytic conformal transformations. 
This includes besides the sum over j G X a summation over the multiplicity n,j a b of these 
states for fixed a,b,j, i.e., we can think of the Vir representation spaces Vj $ as being 
labelled by pairs (j, The "dual vacuum state" is defined by a boundary field placed 
at infinity hm^oo c a ,b,j;p,/3' x 2Aj (0| a ^*,/3'( a; ) where c a ,b,j;p,p' is a normalisation constant 
and (3' = 1, - ■ ■ , rij*b a is an index of type f * fe ) . Accordingly the trace of the operator e~ THab 
computed imposing the periodicity w ~ w + T in time direction can be written as a sum 
over (j,(3) of characters Xj,p > with the summation over j3 leading to fl2.15| ), with a and b 
exchanged. 

Boundary fields b ^j^(x) appear as ordinary fields ^j(x) decorated by a pair of indices 
a, b according to some rules. In the limit L/T — > oo their 1-, 2- and 3-point functions 
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are given by the corresponding invariants of tyj with respect to sl{2, R) with normalisation 
coefficients depending on a, b. As for the ordinary fields the 1-point function is non-zero only 
for the identity operator 

(0\ b *l a (x)\0) = 6 n 6 ba 5 ala (l) a (4.8) 

with the restriction on a, b coming from n\ a b = 5 a b- As for the ordinary CVO the product 
ai d ^'j a 2 °f ^ wo boundary fields is defined only for coinciding c = d. Similarly the initial 
and the final indices in a vacuum expectation value of a product of boundary fields, are 
restricted to coincide (due to the periodicity in the strip time direction, the boundary half- 
line being effectively closed) but in distinction with the ordinary CVO they can be arbitrary 
and not just equal to the identity 1. The 2- and 3-point functions read 

(~<ab 

3/19 J 



/~iabc 

<or*u(*i) b n^) °*i««(^)io>t = s ad ] -^f^f -i , + X2 + x 3 + Xl , 

F12I ij 1^23 1 jk Fsil^w 

(4.10) 

where A& = Aj + Aj — A& and Xij = Xi — Xj. 

The functions ( [4.9| ), ( |4.10[ ) are invariant with respect to SX(2,R), with representations 
denoted by a pair (5, e — ±) , see ^2[ ; here we choose £ = 1 corresponding to taking the 
modulus of the multiplier of the SL(2, R) transformations and the expressions in ( [4.9|j , ( 4.10p 
imply trivial monodromy of the boundary field correlators. In the sl(n)k WZW models the 
fields carry an additional tensor index, or, in a functional realisation, depend on an additional 
(multi) variable X accounting for the representations of the "isospin" sl(n, R) algebra and the 
n-point functions involve also n-point invariants with respect to this algebra. For example, 
in the sl{2) k WZW case the fields tyj(x,X) can be described in terms of a pair of real 
variables p8|, the coefficients in the polynomial expansion with respect to X representing 
the horizontal algebra descendants. In this case the isospin labels are 2j = 0, 1, . . . , k , and 
the 2- and 3-point invariant correlators contain additional factors Xff along with any x a 2Ai . 
For simplicity we adapt the notation for the minimal W n models rather than their WZW 
counterparts, omitting the explicit indication of the "isospin" variables and the corresponding 
invariants. To keep track of the various possible three-point invariants, we shall retain the 
multiplicity index t as in (|4.10|) . 

Up to the normalisation constant and up to phases, (|4.10|) is the 3-point invariant 
function (0|^ 5 x< 4>j,t 0fc,i fc |O) of the ordinary CVO 4>j,t{ x )- Here t is a coupling index of type 
(\ k ) , t = 1, 2, . . . ,Njk l * . Two kinds of permutations act on these couplings, see [Q: 023 

(l?) ~~ > (ji) ' an< ^ ai3 ' (ij) "~ > (ip*) ' ^ or s i m P nc ^y i n the sequel we denote the one 
value indices indicating couplings with one label of type 1 set to 1, like 023 (1*) , l'j , or, l a 
(corresponding to couplings of type Q-J , (j^j , or, (® a J resp.), simply by 1. 



29 



Motivated by the form of the 3-point function, the operator product expansion (OPE) 
of (primary) boundary fields b ^l a (x) is defined according to 
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where P is an index for the descendent states of the representation V p with p £ I. The indices 
cti , ac 2 , j3 account for the multiplicity of vertices of type Q c J , (J a ) , Q a J , respectively, 

a, 6, c £ V, i.e., a\ = l,2...,n ic b , etc, while t is that of a standard vertex (^) , t = 
1,2,... Nij p . We will often restrict for simplicity to the si (2) case, so that the index t can 
be omitted. From the 1-point function 
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(4.12) 



With the normalisation of the CVO 

P i 

the numerical coefficients ^F, 



v 



cp 



i 3 
b a 



implied by the second equality in (|4.11|) 
with i,j,p £ X, a, b, c £ V, represent the OPE 



coefficients of the boundary fields and their determination is part of the problem. They 
are reminiscent of the matrix elements of the fusing (or crossing) matrices F (whence the 
notation here for this "second" fusing matrix), which serve as OPE coefficients of the usual 
CVO |33|. The definition (|4.11|) extends to descendent fields in the l.h.s. as for the usual 
CVO. Symbolically ( [4.11] ) can be written as 



( b ) ( c ) =y {1) f c 



p,p,t 



* 3 
b a 



(3 t 



(4.13) 



and can be depicted similarly as the standard Moore-Seiberg diagrams, see Figure £|. De- 
noting by U h pa the space of boundary fields of type Q b ) we have dim U b a = n pa b while the 
space of standard CVO, Ufp has dimension given by the Verlinde fusion multiplicity Nij p , 
dim Ufj = Nij p . Thus we can interpret as a linear operator 



u 



3 a 



pa ' 



(4.14) 



the dimension of the two sides being identical, according to fl2.22| ). Given the 1- and 2-point 
correlators above the computation of the general boundary field n-point functions is reduced 
to the computation of the conformal blocks of the standard CVO 
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Figure 4: Graphical representation of ( 4.1 3| ) . To stress the presence of two types of vertices, we 
distinguish them explicitly on this particular figure only 



Comparing with the 2- and 3-point functions, see Figure |5](a) and (b), we have 
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Figure 5: (a) and (b): boundary field 2- and 3-point functions 
The 2- and 3-point normalisation coefficients are assumed to satisfy the symmetry 
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conditions 



/~iab /~tba 

u i*i;a/3 ~ ' ii* ;f3a 



s~iabc ribca 

ijk;ai ct2 «3;t jki;cn2 013 ai;ai3(T23{t) 



013 03*) 013 (<**)/ ' 



(4.17) 



i* ;cri3(o^) 013(0!*) <Ti3(af );o-i3(i*), 



The first equalities in ( [4. 17] ) are cyclic symmetry relations, see Figure || while the second 
equalities come from an antilinear ("CPT") transformation, which in particular sends the 
field b ^/ c jf3 (x) to its conjugate c ^ b - ai3 m*) 



' ,,r ' ' —x) with multiplicity indices consistent with nj a c 
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Figure 6: The cyclic symmetry of 3-point functions 
The cyclic symmetry relations imply 
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while the the second equalities in (|4.17l) lead to 
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Combining (|4.18|) and ( |4.19| ), one recovers Q4.16|) . 



(4.19) 



(4.20) 



4.2.2 Bulk fields and bulk-boundary coefficients 

We turn to the second ingredient of the Cardy-Lewellen boundary CFT, the bulk fields. The 
half-plane bulk fields <f>j(z,z) , z = x + iy G H + ,z = x — iy transform under a representation 
of D> H > ]36| realised by differential operators 



H 



L n (Ai,z) + L n (Aj,z) 



(4.21) 
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and characterised by a pair / = (i, i) of weights. In cases when there is more than one field 
with the same labels a more involved notation like / = (i,i;a) is needed, but usually 
omitted for simplicity. For type I theories as well as for arbitrary scalar fields both i,i G £ , 
while in general i, i G Z. 

( H) 

The invariance with respect to the subalgebra spanned by L y ±1 ' determines the 1-point 
function of 3>(i,i)(<z , z) as well as the 2-point function { a ^^ a , e.g., 

(0|X«(*i)*(i.?)CM)l°> ( 4 -22) 

C a 

p,(i,i),a,t . -p. 

x\ > Re z , 



[Z - z)Ai+Aj-A p f Xl _ ^)A,+A P -A. f Xl _ ^Aj+Ap-Ai 



while "^pa) is defined for Kez > x± by the analogous expression with 

replaced by z — xi , z — x%. Requiring the symmetry of this function under the exchange of 
the two fields, i.e., the independence of the ordering, leads to the constraint Aj — G Z. 
The r.h.s. of (|422|) is the 3-point block of the standard CVO (O^^fci) & jt (z) <fc,i ? (2)|0) 



with t a coupling index of type ( ^ ) , t = 1, 2, . . . , N^f* . Consistently with this the (primary) 
bulk field can be represented for small z — z via the decomposition 

^(z ,z) = J2 ^(S) E (P, ^IM* - 0) X«*p(*) (4-23) 

— V a,Q - a vl> a (z) + 

a,a ,p£l ,t V / 

V a ' Q Rf '* - a Vl> a (V) + 

^ T , <*) (2zy)A i +A ? -A J> *p,«W + •- 

which extends to descendents. Here a ^fp a (z) are "unphysical" generalised CVO obtained 
extending to the (full) plane the boundary fields of the previous section. Their OPE is 
determined by the same fusing matrix ^F, i.e., as in ( |4.11| ), the latter extended to complex 
arguments Zi , Rezi 2 > 0, with \xi 2 \ replaced by z 12 ; we shall need only this fusing property. 
The constants a ' a B?.K (" bulk-boundary reflection coefficients") in this decomposition 

depend on two couplings of different types, (j'.jj and Q a J . Note that the coefficients used 
here differ by a phase from the traditionally normalised coefficients [IT|, [0, which will be 
denoted °' Q 'Bfy(C 'L) , i.e., 

°- a B^ = e ^i^+^-^) a > a B^(CL) . (4.24) 
The decomposition ( |4.23|) , symbolically written as 



see Figure 0, reduces the computation of the n-point functions of ^uj) to the computation of 
the blocks of the generalised CVO a ^p a (z), which combined with their OPE (the extension 
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from |a?i2 1 to Z12 of ( |4.11 )) allows to recover all correlators in terms of standard conformal 



blocks. The invariant 1-point function projected onto the boundary state a reads 

($ (iii) (z,*)) a = 5r-i _±% i (l)a. (4.26) 

\Z Z) 

Omitting the trivial indices and simplifying the label (i, i*) to i, one has in particular a B\ = 1 
for any a. 

The OPE of the half-plane bulk fields $( fe ^(2:, z) is defined according to 

= E EOvi^mWI^o) (j,J|0s,f(^2)ir,o> $ {iJ);WJ) (^,^) 

j',i,M J, J 

= E_^%%i)(^^)+-- (4-27) 
The coefficients D^^yfn are related to the full-plane bulk OPE coefficients, see below. 

4.3 Boundary CFT duality relations 
4.3.1 Cardy-Lewellen equations rederived 

We collect in this section the set of equations resulting from the sewing constraints on the 
various OPE expansions ||12|| ; some of these equations can be interpreted as expressing locality 



(symmetry) of the boundary CFT correlators. For simplicity of notation we shall sometimes 
omit the explicit indication of the coupling indices of type Q fc J and the corresponding 
summations, i.e., the equations will be written essentially for the simplest sl{2) case. However 
we shall keep the charge conjugation in the indices of I so that the general formulae can be 
easily recovered. In the si (2) WZW case the braiding phases are given by the shifted scaling 
dimensions Aj Ug — j, instead of A j U9 (since the pair of coordinate and isospin variables is 
moved as a whole). Then formulae work equally well with the same fusing and braiding 
matrices, without additional signs, as for the corresponding subfamily (1, 2j + 1) of fields in 
Virasoro minimal models. 

Applying (|4.11|) in different ways to the 4-point function of boundary fields that is 
demanding associativity, leads to a relation connecting the two types of fusing matrices 
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Figure 7: Graphical representation of the decomposition ( 4.25| ) of bulk fields 
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and F, the fusing matrix for the ordinary CVO, which reads symbolically 

p Wp (i)_p = Wp W_p 



(4.28) 



or, more explicitly, 
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(4.29) 



The identity ( |4.29| ), when restricted to the si (2) case, is a slightly simplified version of the 
equation (L 3.29) in [|T^] and can be also obtained from the latter using the relation (|4.16| ) 
and dropping a (non-zero) factor of type ^F a i. The direct derivation of this pentagon-like 
identity depicted in Figure || is analogous to the derivation of the standard pentagon equation 
for the fusing matrices F since the boundary field n-point blocks are analogs of the ordinary 
(n + 2)-point conformal blocks with an additional constraint due to the delta function in 
the 2-point boundary block. The relation (|4.16 ) is reproduced from the pentagon identity 
(|4.29| ) for particular values of the indices. 

Imposing the symmetry of the 3-point function ( b ^j S (xi) $(j ^(z, z) "'^Jx^)) = 
( b iffj S (xi) a \l/^ 7 (x 2 )$( i j)(z, z) ) one derives following Figure 9, 

1 r • , 1 P' 
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p p 
b b 
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P'P 



j k 
b b 



£ a > a Bfa (l) a UF al 



i ot.at 1 



S S 

a a 



E i7r(2Aj-2A m +A fc +Aj-A p ) p 
8', 



- 7 8 



j i 


F * 


k j 


k* i 


- 1 mp* 


i* i 



(4. 



On Figure ^ the braiding matrices B(±) appear, see Appendix E. In the r.h.s. of ( 4.3U| ) we 
have also used the cyclic symmetry relations ( |4.18|) ,( f4.19|) . Using furthermore these relations 
a factor of type ^F al (and the related summation) can be dropped in both sides of (|4.30|) 
which leads to a slighly simplified version as compared with the original equation (L 3.32) 
in 
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From the 2-point function (<J>o. ^(zi, Zi)<&aj)(z 2 , z 2 )) a using either the OPE formula 
Qgg7| ), or ( ^23|) (follow Figure 10 with i = 1), we obtain 



a,a R r, Sl a ,/3 r>r \s 2 (1) j? 
(k,k) 13 (1,1) ral 



a, 



a a 



la U 



(4.31) 



J a/3 



agl in(A k +Ai-A r -Aj) ^ D (j,J_r>hj; p. 
£-~< j £-~< (k,k)(l,l) l r 



t.t 



k k 
I* I 



- CT230-13(S2) CT23(si) 



J (713(F) t 



4 'Mixed' pentagon identities analogous to (4.28) appear in the framework of 'weak Hopf algebras' as part 
of a 'Big Pentagon identity' j6j|. The counterparts of ^ F are interpreted as kind of "3j - symbols" along 
with the standard interpretation of the fusing matrices F as "6j - symbols" . 
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Figure 8: The "mixed" pentagon identity 



or, equivalently, 



n (j,r)-,t,t a nl „«vr(A fc +Aj--Aj) 

±y (fc,fc)(«,0 -°i e 



E iirAr4,apr,si a,/3 p>r* ,S2 (1) p 



a ,/3 ,r ,si ,«2 





i i' 


' k f 


a a 


a/3 


k* I 



(4.32) 



012 (t) f 
fl3(si) S 2 



Lastly from the 3-point function {$(k,k)( z i) ®(i,i)( z 2) a ^i t -y( x )), we obtain, see Figure [U] 



Ea,a r>r _ a,/3 r>t _ (1) p 



a,/3 



r t 
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(4.33) 
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For i = l (|4.33| ) reduces to ( |4.31|) . The sum over s in the r.h.s. of (|4.33|) represents 
up to phases one of the sides in (an auxiliary) hexagon identity, resulting in permuting 
{l,i,k} to {k,i,l} i.e., can be written as £> 2 3(— )5i 2 (— )£> 23 (— ) and thus can be replaced 
by Biz(— )-B 23 (— )-B 12 (— ). This gives an alternative representation of the r.h.s. of ( |4.33|) 
obtained from the above by replacing everywhere (k,k,j , r) — > (1 ,1 , j ,t*). This is the 
original form of the equation (L 3.35) in []12"| |, when ( |4.24j ) is inserted, with furthermore 
inverse operator ordering convention and opposite overall sign of the phase. 
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Figure 9: Derivation of (KM 



The locality of this 3-point function (the symmetry under the exchange of the fields $) 
implies also 



(k,k){l,T) 



i yk+si-sj n (i,i);o"2_3(t),o-23(F) 



s k := A k - A- k , (4.34) 



while from the associativity of the OPE ( |4.27| ) one obtains in particular the relation 



n U,3);t,t n (i,i) 
^(fc,fc)(M) ^UMrJ*) 



j~ ) (k*,k*);r7 13 a 2 3(t),a 13 (723(t) J") (!.!)_ 

C'»0C3'*J*) ^ {k* ,k*){k,k) ■ 



(4.35) 



All the above equations hold true as well with a sign e = ±1, inserted in the expo- 
nents of all the phases in these equations including ( |4.24|) , and replacing the bulk-boundary 
coefficients B with B e . Thus when rewritten in terms of the Cardy-Lewellen normalised 
coefficients B(CL) all equalities are true for both choices of sign. 

Remark: The Lewellen equations in the diagonal s/(2) case were recently confirmed in [7C]. The 
seemingly different version of equation ( 4.3C| ) in [f4] is in fact equivalent to the original Lewellen 
equation, after taking into account one of the duality relations (a hexagon identity, see ( |E.6[) in 
Appendix E) for the braiding matrices. On the other hand the derivation of the versions of ( [4.31| ), 
( |4,33| ) in |[4[ appears to be affected by a missing phase in the intermediate (and needless) formula 
(17) in [jfj]]. This phase is compensated in the final formulae following from ( |4 . 3 1| ) (like Q4.44| ) 
below) by another phase due to the presumably neglected difference in the normalisation of the 
bulk-boundary coefficients (like in (4.24)) as compared with that in [11, 12]. 
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Figure 10: Derivation of (fO|) 
4.3.2 More pentagon relations 

Before we turn to a discussion on the implications of the Cardy-Lewellen equations we shall 
introduce one more ingredient to the scheme. It is natural to assume that there exists a 
"third fusing matrix", a matrix inverting WF, ( 3 )fWf = I = Wf^F, or m ore explicitly, 
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(4.36) 
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Along with the standard CVOs (jjjj , this matrix involves new "couplings" of type (^ p J , 

which can be thought of as obtained by a permutation cr 23 from the boundary fields (p Q )^ > 

whence the notation (3' = 023 (/3). The matrix satisfies a "mixed" pentagon identity 
analogous to fl4.28| ) 

(3) p (3) p p = (3) p (3)p _ 



6* j 



0-12(03) 012(02) 



(4.37) 



and summing 



it 



Furthermore multiplying both sides of (|4.29 ) with ( 3 'F m / c * 

over c, 02,03, using ( |4.36|) in the r.h.s., we obtain another equation of similar form 

(3) p (1) p (l)p = (l)p p 

which implies various useful relations obtained for particular values of the indices. One of 
them reproduces the inverse property of ( 3 >F, another one reads 



E (3) ^* 



a 3 
a* j 



* "I 0-23(0) 0-23 (P) 



± al 



-111 



J 3 
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J 013(7) "13 03) 



(4.38) 



where dj = Sjx/Su = F u 



dj E (3) ^* 

a,0 



J 3 

3 3 

a j 

a* j 



1 -1 



is the quantum dimension. It furthermore implies 

023(0) 023 (/3) r . -1 1 1 

11 L J 013(0) 0-13 



n 6 



(4.39) 



Using the inverse matrix ^F, as introduced here, the half-plane bulk field can be also written 
as a product of generalised boundary fields (a "bilocal" operator) 



a?a(M)= E ( E ^ (3) ^ 
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a,b,a,f3 p ,7 ,t 



a" 2 
a* i 



" Ol2023(/3) 012023(a) 
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which reproduces the small z — z expansion in ( |4.23| ); compare ( [4.40 ) with the chiral decom- 
position of the (full) plane physical fields 

*5,i,(c,c)= E w ^lfS)^(0®Vf,,(C). 

k,k,l,l,t,t' 

Finally we shall exploit the inverse ^F of the matrix ^F to rewrite ( 4.3U| ) in another 
equivalent form to be used in the next section. Namely we apply the inverse to ^'Fh* in the 
r.h.s. and obtain 
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In the third line we have used the symmetry relation ( |4.16| ). 



4.4 Consequences of the bulk-boundary equations 
4.4.1 The Pasquier algebra and its dual 

In this section we analyse some important consequences of the set of equations derived. We 
start with equation ( |4.41| ) , an inverted version of the first Lewellen bulk-boundary equation 
(|4.30| ), in which we take s = 1 = p. This implies k — j* , i — i*, a — 1 — a' — /3 — (3' . The 
sum over m in the l.h.s. is proportional to the modular matrix Sji, see (|E.9|) , while the sums 



over the coupling indices 5 , 7 are worked out using (|4.39| ) , the final result being 

Sji 



a B] (1>« 
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b 



11 



b b 



.1 " 



B] (l) t 



(4.42) 



For simplicity we have done this computation in the sl(2) case but it extends straightfor- 
wardly to arbitrary rank leading to the same formula. 

Comparing ( }4.42| ) with ( |3.18| ), we see that it can be identified with the realisation ( |3.18| ) 
of the relation ( |3.4| ) in terms of the eigenvalues 7 a (^) = of the graph algebra matrices N a . 
Namely we can identify the ratio h B\ (l) h / a B\ (l) a with the ratio Recalling the 

expression for (l) a in ( |4.(j| ) we find a relation between the boundary state coefficients ip l a and 
the bulk-boundary coefficients a B] 



a B l = e inA, *Bl(CL) 



fa 



I a 



a* 



d: 



(4.43) 

Conversely, if we assume 



where for the time being Cu* is an arbitrary constant, Cn = 1 
the identification ( |4.43| ) (as, e.g., derived in the sl(2) case by other means in |p"T 



with 



1, see also 



we recover the relation 



or, 



directly from one of the 



bulk-boundary equations. As discussed in Section 3.5, from this relation we reconstruct the 
graph algebra. 



In the diagonal case £ = X, where ip J a = S a j, the relation ( |3.18|) coincides with the 



Verlinde formula, i.e., the standard fusion algebra realised by its characters. On the other 
hand the Verlinde formula is known [l4j] to be recovered from the diagonal version of the 



other bulk-boundary equation, the Cardy-Lewellen equation ( |4.31| ) to which we now turn. 
This equation simplifies for r = 1, leading to k = k* , I = I*. Using ( 4.43| ) and denoting 
Pi(a) = ip l a /i/jl the equation ( |4.31| ) turns into 



p k {a)pi(a) = M M J pj(a) 



(4.44) 
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where 
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(4.45) 



<Tl3(t) * 



and Mfc/ J = 
ing O, 



if the corresponding Verlinde multiplicity vanishes. Alternatively, invert- 



l kl 



aGV 



k,l,j e£. 



(4.46) 



Let us first look at the diagonal case in which according to ( f4.46|) the constants M\J 
coincide with the Verlinde fusion rule multiplicities . This is confirmed also directly 
by the alternative expression (|4.45|) provided by equation (|4.31| ) as we shall now show. In 
the diagonal case, denoting D^ k }^^ = Cq. we can use ^ e mver t e d equation (|4.32|) 



taken for a = 1, a choice which trivialises all summations, since n r \ l 



N, 
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result (pointed out in the si (2) case in [[TO]) 
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Taking in particular j = 1 ( 4.471) gives 



r ,2niA h 
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^ ^ (k,k*)(k* ,k) ■ 



(4.48) 



Comparing with ([4.43 ) taken in the diagonal case we see that we can identify the undeter- 
mined constant Cu* with the normalisation constant of the bulk 2-point function. We shall 
retain this identification of Cu* in the non-diagonal cases (at the same level as the given diag- 
onal case) which amounts to setting the relative 2-point normalisation to 1. Combined with 
( 4.35Q and ( 4.43[) the relation ( [4.47D leads to a symmetry of the fusing matrices analogous 
to the cyclic symmetry ( 4.18Q 
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k* 



J l l 
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(4.49) 



kt and completes the argument; alternatively the same conclusion is 



Inserting ( [4.47] ) back into (|4.45|) and using ([4.49 ) reduces the sum over t to the standard 
pentagon identity specialised for some choice of the indices (cf. the analogous relation (|4.38| )). 
Finally we are left with the sum over the coupling index t = t^ kl ^j , which reproduces the 
Verlinde multiplicity Nj 
achieved using the simple choice of gauge 

Note that the relation fl4.47|) , with ( |4.48|) accounted for, is a linear version of the standard 
(quadratic) relation for the full plane diagonal OPE coefficients which results from locality of 
the (full) plane bulk fields 4-point functions, see Appendix E. In the sl(2) case this identifies 
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the OPE coefficients of the half- and full-plane diagonal bulk fields. The identification 
extends to the nondiagonal si (2) scalar OPE coefficients as can be seen generalising to 2- 
point bulk correlators the computation of the limit L/T — > oo of the 1-point correlators in 
llj, leading to ( ggg ). 

In the general (non-diagonal) sl(2) cases characterised by a fixed level (central charge) 
we can express the fusing matrix in the r.h.s. of (|4.45 ) in terms of that in the r.h.s. of (|4.47|) . 
Using once again the si (2) versions of the identities just described, we express it in terms of 
the diagonal OPE coefficients at the same level, obtaining for = 1 

W = <$5xM) = D Uw) l c iAm ■ ( 4 - 5 °) 

The relative scalar OPE coefficients d^f\nn have been computed for the si (2) WZW and 
the Virasoro (unitary) minimal models, see, e.g., for an exhaustive list of references. 
Now using the expression for the eigenvectors fl2.38|) , they can be computed as in ( 4.46Q for 
all minimal models. 

The matrices (M^)^ = can be seen as a matrix realisation of an associative com- 
mutative algebra with identity, distinguished basis and an involution *. In ( 4.44Q the algebra 
is realised by its 1-dimensional representations (characters) given by ratios of elements of 
the eigenvector matrix defining rij. This algebra, traditionally called the "Pasquier algebra" 
("M" -algebra), is dual in the sense of ref. JTIJ to the graph N- algebra considered in Section 
3 but unlike its dual, its structure constants are not in general integral, but rather algebraic 
numbers. In the simplest si (2) case the squares (Mm-*) 2 of these constants are rational num- 
bers for all A-D-E cases; this rule persists for most of the si (3) cases but is broken by two 
of the graphs £^ and £2 corresponding to the exceptional modular invariant at level 
k + 3 = 12, see Appendix D. The type II sl(2) cases, D odd and E 7 are again distinguished by 
the fact that the sign of some of the multiplicities M^f is negative and this is a basis inde- 
pendent statement in the sense that there is no choice of basis to make all M 3 kl non-negative, 
contrary to the Type I cases D even , E 6 ,E 8 , and this is a general feature of type II theories. 

The formula Q4.50D extends beyond the si (2) case for (k,l,j) such that = 1, i.e., in 
cases with trivial Verlinde multiplicity the matrix elements provide the (relative) OPE 
coefficients d9ffk*^ni*y For non-trivial Verlinde multiplicities > 1 the relation between 
the constants ( |4.46|) and the OPE coefficients is not so direct. Let us give a si (3) 

WZW example which illustrates the relation ( |4.45|) . There are three graphs found in [|3£j 
which correspond to the exceptional block-diagonal modular invariant at level k + 3 = 12, 

(12) 

see Appendix D, where these graphs are denoted by E\ , % = 1,2,3. One can pick up 
triplets of weights (i,j,l) such that the Verlinde multiplicity of the diagonal sl(3) model 
at level k + 3 = 12 is trivial, iVy = 1 and check the values of the corresponding Pasquier 
algebra structure constants Mi/ for each of the three graphs. The result is that, comparing 
in particular and £3 , there exist such triplets leading to different values of Mi/ for 
the two graphs. Since for trivial Verlinde multiplicities the formula ( [4.45 ) gives a direct 
relation between the two types of constants, Mi/ = d^l'^'^y this result suggests that 
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there are two different solutions for the bulk OPE coefficients in this case. Only one of 
these two non-diagonal solutions, namely the one which can be associated with the type I 
graph S[ 12 ^ was recovered in |53|, exploiting a set of equations for the M -algebra structure 
constants. This set was derived from the bulk CFT locality equations assuming an additional 
(quadratic) constraint on the OPE coefficients in theories with an extended symmetry; some 



of its consequences were also reproduced in the abstract framework of |55| , in particular the 
relation nn a = mult a (i) discussed in Section 3. Precisely this relation fails (and hence the 
assumptions on the OPE coefficients in ||53|| ) for the graph £3 , which otherwise satisfies all 
the requirements of type I. 



one can 



We conclude with a comment on the OPE coefficients. As discussed in |25|, |26 
relate in the limit L/T —>■ 00 the correlators of the half- and full-plane bulk fields ^f(z, z) 
and $^(C,C); looking at the two dual representations of the partition function with field 
insertions; in particular ( f4.43| ) was recovered in this way. Though this transformation needs 
to be elaborated for higher rank cases it seems reasonable to expect (and in agreement with 
( |4.47| )) that using the two choices of the automorphism Q, discussed in Section 2.1, we can 
identify in this way the OPE coefficients of the two bulk fields with either I' = I = or 
/' (./•./•)• 



A bit of history: The algebra ( |4.44| ) defined through the eigenvectors of the A-D-E Cartan 
matrices first appeared in the context of the sl{2) A-D-E lattice models proposed by Pasquier 
[50] a short time before the Verlinde fusion rule formula (the "A" algebra in the sl{2) case) was 
found. The interpretation in terms of a pair of dual C ('Character')-algebras was proposed in 
|5l| in the discussion of the set of graphs found in |3£| as a generalisation of the Dynkin diagrams 
associated with the modular invariants of sl(3) WZW and minimal models. The fact that the 
relative scalar OPE coefficients <$ff/Ln ^ of all A-D-E series of the sl(2) WZW (or the subfamily of 
fields (1, s) in unitary minimal models) coincide in a suitable basis with the Pasquier algebra (4.44) 
structure constants Mj-i 1 was the main result of \d2\. It was established through a case by case 
check, supported by a lattice model derivation in which the same coefficients appear considering 
representations of the Temperley-Lieb algebra. CFT locality constraints resulting in formulae 
quite similar in spirit to ( j4.45|) were furthermore exploited in [63], [54| as an ingredient in the 
construction of generalised Pasquier algebras and thus of new examples of graphs related to sl(n) 
modular invariants, extending the results in |p9| , [pi]] . The authors of [52], [53|, |34|] were however 



not aware of the parallel development of boundary CFT, and in particular of [11], [12], where the 
equation ( 4.31| ) first appeared. The importance of the algebra obtained from this equation at r = 1 



was recognised and stressed in [14], where a representative example of the sl(2) WZW -D dd series 
was considered, for which the set V and the characters Xi (°0 °f the algebra were explicitly described. 



Presumably the authors of [14] were not aware of the general A-D-E result in [52|. In the same 
framework of boundary CFT the Pasquier algebra reappeared recently in a systematic study of 
orbifold theories, see [f29f | and references therein, under the name "(total) classifying algebra". 
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4.4.2 Relation to the Moore-Seiberg set of duality equations 



We have seen that the two Lewellen bulk-boundary equations ( 4.30Q and (|4.31 ) when re- 
stricted to some particular values of the indices p in b B( i ^ become in some sense "dual" 
to each other, recovering the two dual C-algebras, the graph and Pasquier algebras. These 
algebras are identical in the diagonal case, reproducing the Verlinde fusion algebra, which 
suggests that in this case the above two equations might be related. 

On the other hand let us recall that the original derivation p4| , [33|] of the Verlinde 



formula relies on the use of one of the basic Moore-Seiberg duality relations, namely the 
equation resulting from the modular property of the two-point functions on the torus, see 
(|E.7|) . It involves the fusing/braiding matrices F or B and the modular matrix Sij(p) (in 
general 5|'-*(p)) for the 1 -point functions xf\ r : z ) on the torus, the index p standing for the 



(representation) label of the inserted CVOs (v) / and (^/j , log z' = (log z)/r , see |53 

The alert reader may have already noticed the full analogy between the Moore-Seiberg 
torus identity ( |E.7| ) and the second version ( |4.41| ) of equation (|4.30|) . It suggests that the 
quantity taking over the role of the modular matrix S(p) is the bulk-boundary reflection 
coefficient B p . In the diagonal case this correspondence is precise, i.e., the two are identical 
up to a constant. 

Indeed first note that equation (|4.41|) , still considered in the general (non-diagonal) case, 
simplifies for s — 1, that is k = j* ,i = i*, a' = 1. Inserting in the first line the expression 
for the modular matrix Sij(p), see (|E.8|) , and using in the third line ( 4.38| ), we obtain an 
expression for the modular matrix S(p), 



V 



E "Ufa (l) b E {1) F bj 

b,p s 



i s 



(4.51) 



J/35 



Let us concentrate now on the diagonal case X = S. The sums in the r.h.s. of (|4.51|) 
can be reduced to one term choosing a = 1 and using ( |4.12| ), since riji = Nji h = 5j b , 5 = 1. 
Alternatively, one can take a = 1 directly in the original equation ( |4.3U| ) - the resulting 
(linear) formula for B p in terms of the fusing matrices F (instead of the formula for its 
square derived in was first explicitly written down by Runkel |f70| . More explicitly we 
have in the si (2) case, 



l B} 



di dj 



p p 



$ij(p) 
S 



ii 



Sjijp) 



(4.52) 



Both (|4.51|) and ( 4.52|) easily extend beyond the sl(2) case. In particular restoring all coupling 
indices the latter formula reads in general 



3,0p>P*,t 
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ji 



(4.53) 
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With the help of one of the consequences of the pentagon identity, ( [L53| ) can be inverted 
and brought into a form analogous to that of ( 4.52Q . 

The coincidence of the two seemingly very different quantities, the coefficients B p in 
the expansion ( |4.23|) of the half-plane bulk field and the modular matrix S(p) of the torus 
1 -point blocks, is quite surprising and needs a better understanding. We were led to this 
observation trying to find a connection between the two duality schemes, the one of Moore- 
Seiberg involving the torus, the other, of Cardy-Lewellen, involving the cylinder. Thus to 
bring equation ( |4.41|) , derived from the first of the Lewellen bulk-boundary equations (|4. 3Q| ) . 
into a form identical to the original Moore-Seiberg torus duality relation we furthermore 
need to identify the three fusing matrices, , F , ( 3 'F , i.e., 



bp 



k j 
a c 



7 * 



a f3 



* bp 



k j 
a c 



7 t 



a 13 



pb 



c 3 
a k 



a' /3' 



7' t> 



pb 



c 3 
a k 



a' (3' 
7' t' 



(4.54) 



This identification is consistent since in the diagonal case I = £ both mixed pentagon 
identities (|4.29|) and ( [4.37] ) then become the ordinary Moore-Seiberg pentagon identity; see 
also JHJ, where the identification of the boundary field OPE coefficients with the fusing 
matrices in the si (2) case was first established by a more elaborate argument. 

It is now straightforward to show that the first and second lines in ( j4.41|) reproduce 
the two sides of the Moore-Seiberg identity ( E.7|) . Taking into account ( |4.52j ) and (|4.54j) 



the second Lewellen bulk-boundary equation ( [4.3 1|) is seen also to be a consequence of the 
first, i.e., of the same Moore-Seiberg torus duality relation. To show this, one has to insert 
in ([4.31 ) the expression (|4.47|) for the OPE coefficients and to compare the equation with 
( |4.41| ) with p — 1, see also Appendix E. 

We thus see that in the diagonal case the two basic bulk-boundary equations (|4.30|) , 



( 4.31] ) are not independent and are equivalent to one of the basic Moore-Seiberg duality 
relations. The third bulk-boundary equation (|4.33|), a more general version of (4.31), is an 



identity which involves again only the Moore-Seiberg duality matrices F , B , S, and thus can 
be expected, following the completeness argument of f33[, to be derivable using the basic 
Moore-Seiberg duality relations. This in particular implies that any solution of the set of 
Moore-Seiberg (chiral) duality relations provides a solution of the diagonal boundary CFT 
equations. 



Remark: Rewritten in terms of S(p) the diagonal case Lewellen equation (4.31) can be also 
interpreted as a generalised Verlinde fusion formula with (non-integral) "multiplicities" J- qk given 
by some particular F matrix elements. The matrices J- q are "diagonalised" by S(p) with the usual 
eigenvalues S q i(l)/ Su(l). Because of this they realise another representation of the usual Verlinde 
algebra. This formula, which derives from the Moore-Seiberg torus duality identity, appears to 
have been already considered, following a different motivation, in [ff^| . 

We conclude this section with a few comments on the general non-diagonal cases. The 
Cardy-Lewellen boundary CFT can be looked at as a purely "chiral" alternative of the 
usual CFT approach in which we combine left and right chiral blocks imposing consistency 
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conditions. It has its 'price' in that everything effectively "splits" - the set X is replaced by two 
"dual" sets V and 8 (for type I, while for type II we have to retain the whole X to describe non- 
scalar fields), there are two representations of the Verlinde fusion algebra and a related new 
fusion algebra (at least in type I cases); there are two types of "chiral vertex operators", new 
duality matrices, in particular a second fusing matrix and its inverse, ^F, along with the 
standard F, satisfying new duality relations, the mixed pentagon relations, generalising one 
of the basic genus zero polynomial identities; instead of one relation involving the modular 
matrix S(p), there are two independent relations - the two bulk-boundary equations in which 
the role of S(p) is taken over by the reflection coefficients B p . It remains to find a consistent 
solution of the equations at least in the sl{2) case. Some of the ingredients are already 
known and have been recalled above. In particular the solution for the .D-series has just 



been obtained by Runkel [82 



5 . Conclusions and Outlook 

In this paper we have reexamined various aspects of boundary effects in RCFTs. We have in 
particular analyzed the consistency conditions that determine the allowed boundary states 
and fields and their characteristic data, OPE coefficients, etc. We have seen that boundary 
conditions are naturally associated with a graph, or a collection of graphs, whose spectral 
properties (eigenvalues) encode the diagonal spectrum of the bulk theory. This legitimates 
empirical observations made previously on the role of graphs in the classification of RCFTs. 
We have seen that the torus partition function may be fully reconstructed from the infor- 
mation contained in these graphs. We have seen also that in several cases (sl(2), sl(N)i 
theories), this approach provides a substantially simpler route to the classification of RCFTs 
than the study of bulk properties (modular invariants. . . ). We have finally seen that further 
important information about some boundary effects (^-factors, boundary structure con- 
stants) is also encoded in the spectral properties (eigenvectors) of these graphs. The bottom 
line of this analysis is that a triplet of matrix algebras (rij, N a , Mi) plays a central role in 
the whole discussion. These algebraic structures have been also confirmed by the detailed 
analysis of the basic equations of the boundary field theory. In the diagonal case the triplet 
of algebras reduces to one, the Verlinde fusion algebra. Accordingly, we have seen that in 
this case the basic boundary CFT Lewellen equations can be identified with a set of genus 
(the pentagon) and genus 1 duality identities of Moore-Seiberg. This leads to an identifi- 
cation of some of the basic notions in the two approaches, namely, the boundary fields OPE 
coefficients and the bulk-boundary reflection coefficients B p , with the chiral CFT fusing 
matrix F and the modular matrix S(p), respectively (see the text for precise formulae). 

The more general representations rij of the Verlinde fusion algebra and the dual pair 
{N , M} of associative, commutative (semisimple) algebras have been introduced in earlier 
work on bulk (and later on boundary) conformal field theories, but it seems to us that 
the consistency of the whole scheme now appears in its full generality and that boundary 
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RCFTs reveal these features in a simpler and more compelling way than in the bulk. In a 
loose sense, the boundary effects expose better the underlying chiral structure of the theory 
and its algebraic pattern. This should certainly not come as a surprise, as this is in the same 
spirit as the old connection between open and closed strings. 

The study of a RCFT through its boundary conditions, its algebra triplet, etc, still 
requires a lot of work. The derivation of the Cardy equation relies on a technical assumption 
that has been only partially justified, namely the proper definition of unspecialized characters 
with linear independence and good modular properties for general chiral algebras. Also it 
would be good to have a better understanding of the completeness assumption: given a 
certain number of boundary conditions satisfying the Cardy equation, is it obvious that 
we may always supplement them into a complete set in the sense discussed in Section 2? 
Then many questions have been only partially treated: Justify in full generality the validity 
of expressions ( |3.11| , |3.12|) which have been established so far only for particular cases; 
understand better the nature and fusion rules of "twisted" block representations that appear 
in this discussion; set up a general scheme for the systematic classification of integer valued 
representations of fusion algebras; set up with more rigour the formalism of generalized chiral 
vertex operators, their fusing matrices and the ensuing duality equations as a consistent chiral 
approach, alternative to the Moore-Seiberg scheme, etc, such are some of the outstanding 
problems that are awaiting a proper treatment. 

Also it remains to see how our discussion of boundary conditions must be generalized 
in theories where there is no choice of a common diagonalising matrix ip J a leading to N 
algebra with integer structure constants. In the approach of [^, |5£| , in which the numbers 
N are integers, one has to drop the axiom of commutativity, replacing this algebra by a non- 
commutative structure. In that respect, a better understanding of the relation of our work 
with other more abstract approaches — Ocneanu theory of subfactors, weak Hopf algebras 
— would be most profitable. 

Directions for future work also include the discussion of other cases: rational or irrational 
theories at c = 1, c = 2, or M = 2 superconformal theories are particularly important 
cases in view of their physical applications to condensed matter or to string theory. The 
generalization to other types of twisted boundary conditions along the cycle of the cylinder, 
as examined recently in J83|, might constitute another useful approach. Finally the parallel 
discussion of these boundary conditions and algebraic structures in lattice models should be 
extremely instructive and will be the object of a forthcoming publication. 



Appendix A: The Cardy Equation 

In this Appendix, we rederive the Cardy equation (Section 2.2) in the presence of sources, 
which have the effect of introducing unspecialized characters in the partition function. We 
restrict to a conformal field theory with a current algebra. Let {J a } denote the generators 
in the Cartan subalgebra, and v a be "charges" coupled to them. We consider the theory on 
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the cylinder L x T of Section 2.2, call w = u + iv the local variable, < v < L, u periodic 
of period T, and modify the energy-momentum tensor T(w) into 

T\w) = TM-^£, a J«M-f£(^) 2 (A.l) 

As an elementary calculation shows, the last term is dictated by the requirement that T' 
satisfies the conventional OPE of an energy momentum tensor. The central charge is not 
affected by the additional terms. 

One then computes the evolution operators in the two channels of Section 2.2, see 
Figure |l[ For the cylinder, mapped to the plane by £ = e - 27r ™/ T ; the Hamiltonian reads 

H cyl = — [' T du(T'(w) +T'(w)) = — £dw(T'(w)+T'(w)) (A.3) 
2n Jo 2ix J 

= (^(ip+ir-^+E^jr*-^*)) (a.4) 

Note that the additional term in (|A.l ) and (|A.2|) has not contributed to the integral over 
a closed cycle. Taking into account the fact that on boundary states L = L and 
Jg P ^ Q = — Jg P ' )a , we find that the first expression of the partition function reads 

Z b[a = (b\e-^-^^\a) = (A.5) 
where we have defined 

Xj{q,z) :=tr Vj q Lo -^e 2ni ^ ZaJ ° , (A.6) 

and as above q = e~ 47r ^, hence f = 2i^. 

In the other channel, the time evolution on the strip is described by the Hamiltonian 

Hba = ( / 7T-T>(w) + / -—T {w)j (A.7) 

K Jo Alii JO 2711 ' 

and upon mapping on the upper half plane H by z = e 71 "^ we find 

H ba = -(4^--)-— E^ Q --^E(— ) 2 (A.8) 



where now the additional piece in ( |A.1|) contributes the last term. Since the theory with en- 
ergy momentum tensor ( A.l]) and QA.2 ) has still the same operator content and multiplicities 

n i*a = n ib a as before, we may write 

Z bla = tre~ TH »« = e 2 ^^l ^) (A.9) 
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with q = e w l . We then use the modular transformation of unspecialized characters (see 
page 264): 

Xi (q, v) = e 2 ^ § £ vf) (A.10) 



together with the linear independence of the Xi(?> z ) to conclude that ( |2.16|) is indeed true 
in full generality. 



Appendix B: A-D-E Diagrams and Intertwiners 

In this Appendix we establish notations on A-D-E Dynkin diagrams and on the associated 
intertwiners. 

Let G be a Dynkin diagram of the A-D-E type with Coxeter number g. It has n nodes 
that may be coloured with two colours, i.e. its n x n adjacency matrix G a t connects only 
nodes of different colours. This matrix is symmetric and it may thus be diagonalized in an 
orthonormal basis. We call this orthonormal basis ip™, it is labelled by the node a and the 
exponent m (see Figure 2 and Table 1). Hence 

£G V 6 s ' = 2cos — . (B.l) 

b 9 

The ^'s satisfy orthonormality conditions, namely 

E^Vf* = <W S ', S "GExp(G) (B.2) 

a 

E <4>i m = Sat (B.3) 
s'eExp 

Because of the 2-colourability of G, one may attach a Z2 grading r to each node a. One 
proves that if s' is an exponent, so is cr(s') = g — s' and the may be chosen to satisfy 

r js') = ( _i r M C ' . (B.4) 

Moreover all graphs having even exponents, v.i.z. the A, D ^ and E$ diagrams, have 
a Z 2 automorphism 7 acting on their nodes and preserving their adjacency matrix (i.e. 
Ga = G 7 ( a ) 7( - fc - ) , this is the natural Z2 symmetry of these graphs) such that 

<a) = (-l) T(S 'Vf ■ (B.5) 

Finally, one may find in the graph G a distinguished node labelled a = 1 such that ip™ > 
for all m. This special node is typically an extremal node, i.e. the end of a branch for the 
A-D-E graphs; this is generally the end of a long leg, but for the -D dd graphs, for which we 
must choose 1 as the end point of one of the two short legs. 

We list hereafter the explicit expressions of eigenvectors of the various Dynkin diagrams. 
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£ +1 series are the simplest examples of orbifold models. Their fundamental graphs 
can be obtained by folding the A g _i Dynkin diagram so that the nodes aj = a g _i G Db +1 , i = 
1, 2, 3, . . . , | — 1 are identified with the orbit {i} of i under the Z2 automorphism a, a{i) = g — 
i, while the fixed point i = g/2 is resolved into two points az t ± on the graph, denoted n,n — l 
in Figure 2. This implies for the adjacency matrix elements G aia . = A^j + A ia ^ = Aij , for 
i, j 7^ I , and G ag ±a . = Agj and allows us to determine the eigenvectors ipl of G in terms of 
the eigenvectors S^- of the A adjacency matrix. To simplify notation we shall use sometimes 



Eigenvectors of the D21 adjacency matrix 



ijj' 



3 

2 



V2Sij, 

V2 ' 3T 2' 



2 > e 



(B.6) 



fe + ee'W(-l)<). 



For z -odd the orbits {i} belong to T e:c * and can be identified with the subset T = {a^ , i 



1,3, 



' 2 



— 2 , ± } . The matrix Sffy^y is the extended theory modular matrix. The ex- 



pressions ( |B.6| ) can be rewritten in the compact form 

= Sa{j} 



s 



1/ 



\ Qext 



(B.7) 



where S a {jj is a rectangular matrix coinciding for a e T with •S'afj}' while for a = T, 
Soi{{l,±)} = % = and S ai{j} = E; e « = 2^ for j ^ #/2. 



Eigenvectors of the D21+1 adjacency matrix 

Vi = (-l^v 7 ^-, 4=0, z^|; (B.8) 

^§,± ^ ' ^ 21 Jg" JT 2' ^2 

The identity node is chosen to coincide with one of the 'fork' nodes 1 = as. ,+ (denoted 

by L in Figure 2) so that the dual Perron-Frobenius eigenvector tp{ — ipi , has positive 

2 

9 

entries (while ipii = 0). The 'fundamental' node / is identified with ^ — 1, i.e., G = Nz-i. 
Also a* = a for all a, while 7(aa j+ ) = 

Next we display the eigenvectors of the exceptional E r Dynkin diagrams as a matrix 
{ipi}, with the row index a running over the nodes, following the numbering of Figure 2, 
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and the column index j over the exponents in the same order as in Table 1. There too, SV,- 
denote the eigenvectors of the diagonal graph adjacency matrix A with the same Coxeter 
number. 



Eigenvectors of the Eq adjacency matrix 
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are determined from ^ = 



Eigenvectors of the Ej adjacency matrix 
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where a ,b ,c ,d ,e , f are determined from ip{ = ^ S\j Y,ie P Sij where p = {1,9,17}, and 
■02 = ff 1 (The values in the 5'th row come from ipl = V2Sqj for j = 1, 5, 7 .) Explicitly, 

a = [18 + 12^ cos ^]-3, d = [12(1 + cos f 

6= [18 + 12v^cos^f]-5 e = [12(1 + cos (B.ll) 

c = [18 + 12^ cos / = [12(1 + cos ^)]"3 . 
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Eigenvectors of the E% adjacency matrix 
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(B.12) 



where a ,b ,c ,d ,e , f , g ,h are determined from ip{ 
j G S , p = {1 , 11 , 19 , 29} . Explicitly 
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(B.13) 



75 - 30^ 



To such a graph G, one then attaches matrices Vi as follows. The case of reference is 
the Ag-i diagram of same Coxeter number g as G. For this A graph, both the nodes and 
the exponents take all integer values in {1, • • -g — 1}. The if)'s are then nothing other than 
the entries of the (symmetric, unitary) matrix S of modular transformations of characters 
of the affine algebra at level g — 2 
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— sm 
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(B.14) 



in terms of which the fusion coefficients N^ 13 may be expressed through Verlinde formula. 
Note also that 



N . . 9-1.2 

1 y g—ni 



N- 12 



(B.15) 



because of property ( |B.4j ) applied to = S. 

We now return to the graph G of Coxeter number g. The fused adjacency matrices Vi 
with % — 1, . . . , g — 1 are n x n matrices defined recursively by the s/(2) fusion algebra 

Vi = V 2 V-! -V-2 2<i<g (B.16) 

and subject to the initial conditions V% = I and V 2 = G. (One may see that V g = 0). The 
matrices V, are symmetric and mutually commuting with entries given by a Verlinde-type 
formula 
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(B.17) 
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Regarded as (g — 1) x n rectangular matrices, for a fixed, the V, a b intertwine the A and G 
adjacency matrices 

Y,A?V ila b = Y, V ia G v b ■ (B.18) 

V b' 

Regarded asnxn matrices, the Vi satisfy not only their defining relation ( P-16Q but also 
the whole si (2) fusion algebra ( |2.22| ) 



V h V i2 =J2 N hn i3V is ■ (B.19) 



13 



From their recursive definition and initial conditions, it follows that the entries of the V 
matrices are integers. What is not obvious is that these entries are non-negative integers. 
This follows either from a direct inspection or from an elegant group theoretic argument due 
to Dorey |74j] . We refer the reader to [[JpJ, JP| for the explicit expressions of these intertwiners. 



As a consequence of the existence of the automorphism 7 defined above ( |B.5[ ), a Z2 symmetry 
on A, D odd , E 6 graphs and the identity for D even , E 7 ,E 8 , one has 

V- 7(6) = Vj . (B.20) 

Using ( p,16|) and (|3.4j) (i.e. in the notations of this Appendix, V{N a = Y,bVia b Nb) 
one can express the graph algebra N a matrices for all but the -D ev en cases as polynomials 
of Vs with integer coefficients This in particular ensures that they have integer matrix 
elements. For D even one of the extended fusion algebra generators Nff*^ has to be added 
since Va = N e { f +} + Nff_ } while N { = V = V g ^ for i = l,2,...,f - 1. For the three 
exceptional cases E r we have Ni = Vi , i = 1, 2, ...r — 3 , A^ r _i = V r -\ — V r -z , 2 — K- - K--4 > 
N r = V r -2 + V r -4 — V r , which translates into relations between the eigenvalues and given ip{ 
allows to express any %p a in terms of the modular matrix S elements. 

Over recent years, these matrices have made repeated appearances in a variety of prob- 
lems. Originally introduced in the discussion of local height probabilities in lattice mod- 



els 150] and of boundary partition functions B |3^] (see below), they have also appeared in 



the following contexts: 
(i). The "cells" or intertwiners of Boltzmann weights of height models [33, iS| 



(ii). The decomposition of the representation of the Temperley-Lieb algebra on the space 
of paths from a to b on graph G onto the irreducible ones on the paths from 1 to s on 
graph A g _ x |9j according to R a G)b = ® s V sa b R[ A)s . 



(iii). The counting of "essential paths" on graphs [f73][ ; see also recent mathematical work 
by Xu, Bockenhauer and Evans |55|, |56] . 



(iv). The expression of the blocks of the partition function ( |1.2| ) as (p.llj , |3.12j ), see Section 
3.3. 
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The sl(2) intertwiners appear in the computation of the multiplicities m s , of an ir- 
reducible representation b of the finite group, associated with G in the McKay cor- 
respondence | 75| ) 76) , in the SU(2) representations of dimension s' |84l . Namely the 
coefficients of the Kostant polynomials in the generating function Fj, of these multi- 
plicities are given for a non trivial b by J2 C Go c V^ c , where G a b is the adjacency matrix 
of the affine Dynkin diagram and a = is the afline node deleted in passing from the 
affine Dynkin diagram to the ordinary one. The proof of this fact is reduced to the 
recursive relation (IB. 161). 



(vi). These same entries seem to appear ubiquituously in the description of S- matrices of 
affine Toda theories f7^J and in the description of the excitation spectrum of integrable 



lattice models [78, 79, 80 . 



Appendix C: Uniqueness of the Boundary Conditions of 
Minimal Models 

C.l. Matrices with spectrum 7 < 2 

We first recall general results on symmetric matrices with non-negative integer entries and 
with eigenvalues between —2 and 2. 

It is a standard result that symmetric matrices with non-negative integer entries and 
eigenvalues 7 g] — 2, 2[ may be classified. A lemma of Kronecker asserts that the eigenvalues 
are of the form 2cos^ for integers Pi and hi and for the largest one(s), p = 1. One may 
regard any such matrix as the adjacency matrix of a graph. Irreducible matrices correspond 
to connected graphs, and by an abuse of language one may call a matrix bicolourable if the 
graph has that property. One proves that any irreducible bicolourable symmetric matrix 
with spectrum in ] —2, 2 [ is the adjacency matrix of one of the simply laced Dynkin diagrams 
of type A-D-E. 

If one relaxes the assumption of bicolourability, with any symmetric non-bicolourable 
irreducible matrix G one may associate a bicolourable symmetric matrix with a block form 

G' = ( J . The corresponding graph is irreducible and has a Z 2 symmetry that ex- 

\G J 

changes the two colours. Any eigenvalue 7 of G gives rise to two eigenvalues ±7 for G' and 
one thus concludes that G' is of A-D-E type, and its irreducibility forces G' = A 2p . Its Z 2 
quotient G is what we call the tadpole graph T p = v4 2p /Z 2 . 

Finally if one relaxes the assumption of irreducibility, one concludes that any matrix 
(with non-negative entries and spectrum between —2 and 2) is the direct sum of A-D-E or 
tadpole graphs 

G = ®Gi Gi of A-D-E or tadpole type 
and this decomposition is unique, up to the permutation of factors. The uniqueness may 
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be easily proved by induction on the number of terms or on the dimension of the matrix: 
Given a matrix G, one first identifies its largest eigenvalue, of the form 71 = 2cos^-. By 
the previous statement, there is an A-D-E or tadpole graph Gi with Coxeter number h\ 
and exponents rrii, such that all its eigenvalues 2cos^- appear in the spectrum of G. Thus 
G = G± © G", and one may apply on G" the induction hypothesis. If 71 has multiplicity 1, 
this suffices to establish the uniqueness of the decomposition (up to permutations), while the 
case where 71 has non-trivial multiplicity is also easily dealt with. The uniqueness of this 
decomposition implies a property used several times in the text, namely that the spectrum 
(between —2 and 2) determines the form of the matrix up to a permutation of its rows and 
columns. 



C.2. Representatives of n\2 and n<2i 



We now return to minimal models. 
Explicit form of n\ 2 and n 2 i 

It is convenient to work in a basis different from that used in (|2.38 ). In the basis r% = 1, • • • ,p, 
a G G, the second term in (|2.33|) does not contribute to n 12 since N 2pri r2 = 5 ri)2p +i_ r2 = 
for 1 < ri, r 2 < p. Thus n\ 2 = I p <E> V 2 = I p <E> G, where I„ is the p dimensional unit matrix 



n X2 



G 



(G 



V 



G 



\ 



G) 



(CI) 



As for n 2 \, in the same basis, the second term of ( |2.33| ) receives a contribution only from r\ = 
r 2 = p, namely (A^ 2p _i)p P = 1, while V g -\ = T, the matrix that realizes the automorphism 7: 



(C.2) 



Thus one finds that 



n 2 i 



( 

In 

\ 



In 





1 n 

r 



/ 



After conjugation by a block-diagonal matrix with V and I n in alternating positions, which 
leaves the form ( |C.1| ) of ^12 unchanged, n 2 i may be recast in the form 

/or \ 
r r 
■•. ■•. r 
r v ) 

in terms of the tadpole T p adjacency matrix. All the other obtained as universal 

polynomials of the two matrices ni 2 and n 2 i. 



n 2 \ 



\ 



T 



(C.3) 
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Uniqueness of the form of nxi and 7721 

Conversely, suppose we only know that the representation n rs has a spectrum specified by 
the set of exponents £. We want to prove that there exists a basis in which nu and n 2 i take 
the forms and ([031) . 



We first make use of the property that the set Exp(G) is stable modulo the Coxeter 
number g of G under multiplication by any integer coprime to g and is also stable under the 
reflection s — > g — s. We then find that the spectrum of n\2 is made of p copies of Exp(G). 
As explained above in C.I, this implies that in some basis 

rii2 = I P ® G . (C.4) 

For the other generator 77,21, one observes first that the set of numbers that appear 
in flDflP , namely {2cos^}, r' = 1, 3 • • • , 2p - 1, is simply the set {(-l)f +1 2 cos 
r" — 1, 3, • • • , 2p — 1, which is (— times the spectrum of the tadpole T p . According to 
Q2.3H) , this has to be multiplied by (— as s' runs over the exponents of G. Thus if g 
is even (which is the general case except when G = A21) the spectrum of 7721 is made of as 
many copies of that of T p (resp. — T p ) as there are odd (resp. even) exponents in G. For 
A 2 i which has as many even as odd exponents, the same conclusion is still correct! Finally 
one notices that these signs are just the eigenvalues of the V matrix, and one thus concludes 
that 

7721 ~ t p ® r , (C.5) 

where the sign ~ means that it holds in some basis obtained from that of ( |C.4| ) by a simul- 
taneous permutation of rows and columns. From this expression, one can see that 

77 2 i has no row or column with more than two Is, (C.6) 
77 2 i has exactly 77 rows and columns with one 1, (C.7) 

properties invariant under permutations of rows and columns. 

We also know that 77 2 i must commute with T7i 2 . In a basis in which T7i 2 takes the form 
( |C.4j) , 77 2 i may thus be regarded as made of 77 x n blocks that commute with G. We shall 
combine these facts about 7721 as follows: 

• A non- vanishing matrix X with elements in N which commutes with G cannot have 
a row or a column of zeros. Proof: let ip 1 be the Perron- Frobenius eigenvector of G, 
GX = XG implies that Xip 1 is an eigenvector of G with the same eigenvalue, hence 
proportional to ip 1 , Xip 1 = cip 1 , with c ^ since the entries of both X and of ip l 
are non-negative. If X had a vanishing row, Xip 1 would have a vanishing component, 
which is impossible for the Perron- Frobenius eigenvector. If X has a vanishing column, 
one repeats the argument with X T . 

• Any matrix X with elements in N which commutes with G and which appears in the 
block decomposition of r7 2 i cannot have more than one 1 per row or column. Proof: If 
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a matrix X with more than one 1 in a row (respectively column) was a block of 7121, 
because of the property QC.6Q , all the other blocks on the left or the right (respectively 
above or below) of X would have to have at least one vanishing row (respectively 
column), which is impossible by (i) above, or to vanish altogether. In the latter case, 
after a possible reshuffling of rows and columns leaving (|C.4|) invariant, one would have 



either n 2 i 



fx 



Vo 



/ 



or n 2 i 



/ X 

X T 



V '■ '■ 



\ 



(C. 



which would lead to a pair nx 2 , n 2 \ reducible in the same basis. 

It follows that the matrices that may appear as blocks in the decomposition of n 2 % 
must be matrices with one 1 on each row and column, i.e. permutation matrices 
that commute with G. These permutation matrices are the symmetries of the Dynkin 
diagram, and thus are readily listed: 



X 



L r 



I, Ti,i = 1, 



/, r 
I 



if G = A n , D 2q+ i, E 6 
5 if G = D A 

\{G = D 2q ,q>2 
if G = E 7 ,E 8 



(C.9) 
(CIO) 
(C.ll) 
(C.12) 



Here Tj denote the 5 non-trivial permutations of the nodes of the -D4 diagrams, and 
the matrix V exchanges the two end points of D 2q , q > 2. 

One then demands that the symmetric matrix n 2 \ made of such blocks is irreducible 
and satisfies (|C4|) -( |C77|) . This implies that at most one non- vanishing block appears on the 
diagonal. Consistency with the form n 21 ~ T p ® T leaves as the only possibility 



n 2 i 



( 



V 





xl 





x 2 



X, 



X 



p-1 



p-1 
r 



(C.13) 



where Xi, X 2 , . . . , X p _i are chosen among the symmetry matrices of G. A final permutation 
of rows and columns by a block diagonal matrix diag(Yi, Y 2 , ■ ■ ■ , Y p ^i, I) brings n 2 \ into 
the form n 2 i = T p <g> T while leaving the form ( |C.4D of rii 2 unchanged, provided Yj = 
TYj + iXj, hence Yj = r p_: 'Xj_ 1 X^l 2 • • ■ Xj . Then both n X2 and n 2 i have their canonical 
forms Q,0- Q-E.D 
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Remark 

Although it is not required for the present analysis, it may be interesting to look at the 
commutant of matrices of A-D-E type. 

For any G of A-D-E type, with the exception of D even , all eigenvalues are distinct. It 
follows that any matrix X that commutes with G may be diagonalized in the same basis 
as G and consequently be written as a polynomial of G, i.e. as a linear combination of 
/, G, G 2 , ■ • • , G n ~ l . In order to look at cases where entries of X are requested to take values 
or 1 only, and with constraints on the number of l's, it is advantageous to use rather the 
basis of fused graph matrices: X is a linear combination of the linearly independent matrices 
Ni = I, N 2 = G, ■ ■ ■ , N n . The D even case is slightly more involved, since the matrices that 
appear naturally are not independent. 

The commutant of an A-D-E matrix is: 

• A linear combination of the graph fusion matrices N a for G = A, -D dd, E e , E?, E s . 

• A linear combination of the N a and of two of the three matrices H a b, a ^ 6 = 1,3,4, 
that exchange two of the three extremal points of the -D4 graph. 

• A linear combination of the N a and of the two matrices X = T'N 2q and Y = N 2q T', 
where the matrix V exchanges the two end points, for G = D 2q , q > 2. 



Appendix D: s/(3) Modular Invariants and Graphs 

The WZW si(3) theories may be discussed along the same lines as in Sections 2 and 3. Solu- 
tions rii to the Cardy equation are associated with graphs, with specific spectral properties: 
their eigenvalues are given by ratios of elements of the modular S matrix labelled by weights 
of the diagonal spectrum S of the bulk theory. Conversely, spectral properties and the fact 
the n's form a representation of the fusion algebra are not restrictive enough to yield the 
list of possible bulk spectra, as occurred in si (2) (up to the unwanted "tadpole" graphs). 
There are indeed many solutions, i.e. graphs and representations of the fusion algebra, that 
must be discarded as not corresponding to a modular invariant partition function in the list 
(Table 2) of Gannon ||81|| : see |39| for such extra solutions. We may summarise the salient 



features of the analysis as follows, see also the accompanying Tables 2 and 3 and Figures |TT 
and |TJ. 

• At least one graph (or rather one set of n matrices) has been identified for each bulk 
theory, i.e. each modular invariant. But it is not known if this list of graphs and n's 
is exhaustive. 

• Note that the hypothesis of 3-colourability of the graphs that looked natural on the 
basis of the sl{2) case has to be abandonned if we want to cover all cases. This is 
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manifest on Table 3 where it appears that in some cases (namely A^ n '* , and D^ n \ n 
not a multiple of 3, and E^*, the set £ is not invariant under the automorphism o of 
(|C.1J) ) as it should be if the graph was 3-colourable. 

• There are a few pairs or even triplets of isospectral graphs, i.e. different sets of n's 
that give distinct solutions of the Cardy equation for a given bulk theory. These 
graphs/representations should not only describe different sets of complete orthonormal 
sets of boundary conditions for that bulk theory, but also presumably be associated 
with different operator algebras and lattice realisations. 

• In Table 3, which summarises the state of the art, we have also indicated if the graph 
is of type I or type II, following the discussion of Section 3. Some hybrid cases are also 
encountered, in which the M and N structure constants are both non-negative, but 
the N algebra has no subalgebra isomorphic to some extended fusion algebra. 



Notations and footnotes for Tables 2 and 3 

(Shifted) weights of SU(3) A = (Ai, A 2 ) := AiAi + A 2 A 2 , where A 1; A 2 are the fundamental 
weights of SU(3) , A* = (A 2 , Ai) , triality r(A) := (A x - 1) + 2(A 2 - 1) = Ai - A 2 mod 3 . 
Q is the set of weights of triality zero. 
Weyl alcove of shifted level, or "altitude", n := k + 3 , 

Pj5 = {A = AiAi + A 2 A 2 | Ai, A 2 > 1, Ai + A 2 < n - 1} . 
Automorphism a of 

<7(A lj A 2 ):=(n-Ai-A 2 ,Ai). (C.l) 

(a) . One of the two connected parts of the fused Dynkin diagram of type A n _i (the one that 

possesses the exponent (1, 1)). It looks different depending on whether n is even or 
odd. The M, N algebras of A {n) * are positive, as they follow simply from the Verlinde 
fusion algebra N of sl(2). If n is odd, A^* is the connected component of the graph of 
adjacency matrix A 2 n _ x — l made of the nodes a odd (integer spin in si (2)). For a triplet 
of "exponents" A = (I, I), p — (m, m) and p = (r, r) of A {n) % M A ,/ = N lm r + N lm n ~ r , 
and N a b c is the restriction of the Verlinde A n _i algebra to odd a, b, c. For even n, the 
M and N algebras of A {n) * = I + coincide with the Verlinde algebra of Aa-i. 

(b) . The orbifold of A {n \ see |§. 

(c) . The ordinary Z3 fold of AS n \ 

(d) . The unfolded (and 3-colourable) version of A^*. Their adjacency matrix is a tensor 

product by the permutation matrix — ( a 1 ] ; their M and N algebras are simply 

V 1 °J 

obtained from those of A^*, thus also > 

(e) . The Z 3 fold ofS™. 
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Table 2. List of sl(S)k modular invariants; n = k + 3. 



(£)(»)) 



,(") 
++ 

I E AeQnP ^ + ) Ixa + X<tA + X<t 2 a| 2 if 3 divides n 

E AeQnP |" + ) Ixa| 2 + E Ae p|" + )\ Q XaX^- a if 3 does not divide n 

\ J2 XeQnP M ( ELo X^ a) ( ELo X^ a* ) ^ 3 divides n 
E AeQn p^ XaXa* + E AGP (")\ Q XaX*-» tA . if 3 does not divide n 



2 O 
CD 



(E (8) ) Z = |X(1,1) + X(3,3) I 2 + |X(3,2) + X(l,6) I 2 + |X(2,3) + X(6,l) ? + |X(4,1) + X(l,4) ? + |X(1,3) + X(4,3) ? + |X(3,1) + X(3,4) I 

(£( 8 )*) Z = |X(1,1) +X(3,3)| 2 + (X(3,2) +X(1,6))(X(2,3) + X(0,1))* + «■ + |X(4,1) +X(1,4)| 2 + (X(l,3) + X(4,3) ) (X(3,l) + X(3,4) ) * + C.C. 

(E {12) ) Z = |X(1,1) +X(10,1) +X(1,10) +X(5,5) +X(5,2) +X(2,5)| 2 + 2|X(3,3) + X(3,6) + X(6,3) ? 

( E MS) Z = IX(1,1) + X(10,l) + X(l,10)| 2 + |X(3,3) + X(3,6) + X(6,3)| 2 + |X(5,5) + X(5,2) + X(2,5)| 2 + |X(4,7) + X(7,l) + X(l,4)| 2 

+ |X(7,4) + X(l,7) +X(4,1)| 2 + 2 IX(4,4)| 2 + (X(2,2) + X(8,2) + X(2,8) ) X(4,4) + C - C - 
( E MS*) Z = IX(1,1) + X(10,l) + X(l,10) I 2 + |X(3,3) + X(3,6) + X(6,3) P + |X(5,5) + X(5,2) + X(2,5) P + 2|X(4,4) P 

+ (X(4,7) + X(7,l) + X(l,4)) (X(7,4) + X(l,7) + X(4,l))* + C - C - + (X(2,2) + X(8,2) + X(2,8) ) X(4,4) + CX - 
(_E (24) ) Z = |X(1,1) + X(22,l) + X(l,22) + X(5,5) + X(5,14) + X(14,5) + X(11,U) + X(ll,2) + X(2,ll) + X(7,7) + X(7,10) + X(10,7)| 2 

+ |X(7,1) + X(16,7) + X(l,16) + X(l,7) + X(7,16) + X(16,l) + X(5,8) + X(ll,5) + X(8,H) + X(8,5) + X(5,H) + X(ll,8) P 



Table 3. List of modular invariants and associated known graphs for sl(S) 
Footnotes (a)-(e) are explained at the end of Appendix D. 



Modular Invariant 


Graphs 


Altitude 


.Exponents 


Type 


A (n) 


A (n) 


n 


p(«) 


I 






n 


{(i,i)} i<i<L^J 


M,N > 


£>(«) 


V (n) (6) 


n = mod 3 
n / mod 3 


HQ, with (§,f) triple 

/>'".' rg 


M,N > 


£)(«)* 


P W -(Tl23®(4-l*4-l)c (d) 


n 


{(j,j),(n-2j,j),(j,n-2j)} 

i < i < L^J 


M,N > 






8 


(1,1), (6,1), (1,6), (3,3), (3,2), (2,3), 
(4,1),(3,4),(1,3),(1,4)(4,3),(3,1) 




E (s)* 


£(8)* (e) 


8 


(1,1), (3,3), (4,1), (1,4) 


M,N > 


£(12) 


f{ 12) ,* = 1,2,3 


12 


(1,1), (10,1), (1,10), (5,5), (5,2), (2,5), 
and twice (3,3), (3,6), (6,3) 


cf 2) , i = l,3: / 
4 12) // 






12 


(1 1) (10 1) (1 10) (5 5) (5 2) (2 5) 
(3, 3), (3, 6), (6, 3) and twice (4, 4) 


// 


F (12) 


£•(12) 


12 


(1,1), (10,1), (1,10), (5,5), (5,2), (2,5), 
(4,1), (7, 4), (1,7), (1,4), (7,1), (4, 7), 
(3, 3), (3, 6), (6, 3) and twice (4, 4) 


JJ 


e( 24 ) 


£(24) 


24 


(1,1), (22,1), (1,22), (5, 5), (14, 5), (5, 14), 
(7, 7), (10, 7), (7, 10), (11, 11), (11, 2), (2, 11), 
(7,1), (16, 7), (1,16), (1,7), (16,1), (7, 16), 
(5,8), (11,5), (8, 11), (8,5), (11,8), (5, 11) 





-0 vr^) yrQ-Q r.^^ ^Q-Q-C 




Figure 11: The known graphs in the case of sl(3). Conventions : (a) For the 3-colourable graphs, 
the triality r of nodes is indicated by the colour: black r = 0, grey r = 1, or white r = 2; the graph 
represents the matrix n2i if edges are oriented from black to grey, or grey to white, etc. (b) For the 
non-3-colourable graphs, either the orientation of all edges (of matrix 7121, say) is indicated, (T>( n \ 
3 /n, series, £( 8 >), or all links are unoriented (A^* series). 
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24) 



Figure 12: The known graphs in the case of sl(3), (cont'd). 
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Appendix E: Formulae for fusing, braiding and modular 
matrices 



We collect in this appendix some standard formulae for the genus and 1 duality matrices. 
The fusing matrices satisfy several identities implied by the pentagon identity - they can be 
recovered from formulae ( |4.36|) , (|4.38|) , ( }4.16|) in the text making the identification (|4.54|) . 
Symmetries: 



pq 



p q 



k 



pq* 



i* I 
k* J 



p*q* 



I 



Choice of gauge: 



? (o) 
91 



i i* 




d q 


J 3 _ 


a/3 \ 





s 



(E.1) 



(E.2) 



In the sl{2) case denote by the normalisation of the CVO in this gauge. Then for the 

fusion matrix corresponding to CVO normalised to 1, one has 



pq 



k j 
1 I 



\ ^kp^ji 



pq 



^kj^iH^qq* rp(Q) 

\ CP* CP c\ PI 
\ ^ki^jl^pp* 



k j 



or, 



(~1P* (~lP Tp 

^ki^ji^pp* r pq 



^kj^i*l^qq* 1 IP 



k* 
J 



I 

1* 



(E.3) 



(E.4) 



This equation coincides with the quadratic relation resulting from locality of the physical 4- 



point function in the diagonal case. Hence the constants C l kp = C f 



'(k,k)(p,p) 



^kp ^ii* 



^ki* ^pp* 



can be identified with the physical OPE structure constants in this case. For the minimal 
models these constants were computed in flBUfl ; the matrices in the gauge ( |E.2| ) coincide 
up to signs with a product of standard q- 6j symbols, see |[SJ for the latter. 
Braiding matrices: 



B 



pq 



0-23(72) 6 



|9l <T23(P2) 



7rje(A fc +A ( -A p -A 9 ) tj 



pq 



k j 



-| 72 s 



(E.5) 



J/3i/3 2 



The q - analogs of the Racah identity (hexagon identities), e = ±1: 



11 iq 



— 7rieA 



9 F 



qp 



3 7rie(A m +A p — A;— Aj— Aj— A*,) tj 



nip 



(E.6) 



Recall the Moore-Seiberg torus duality identity resulting from a relation in the modular 
group of the torus with two field insertions, namely, S(ji,j2)a = frSfji, j'2) where S^ji, j'2) 
is the modular matrix of two-point blocks, expressed in terms of F and S(p), and a, b are 
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the monodromy transformations moving one of the CVO around the a, b cycles, |33 
Sh(s) E e 2 -( A - A ™) F s * m 

m 

= E SM e^ A »- A n~ A ^ F sg , 



32 * 

fx i*. 

r* ji 

r* j 2 



mp 



3l J2 

i i 



32 Ji 
q q 



Choose s = 1 = r. This implies that j\ = j\ = j* and q = j, hence 





F\ p 


'j 3*' 


J j . 






Fi p 


'j 3* 


J j . 




Sij 




i* i 


F p i 




i* i* 
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The second equality is obtained reversing the sums in ( E.7|) and solving for S r i(s) as above, 
while the third is obtained from the transposed version of (E.7|) taking into account S(p) 2 = 
C e~ mA p. For p = 1 the formula reproduces the ordinary S = S(l) matrix 
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